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We present a study of parity (P) violating contributions to the eigenenergies of stationary 
. systems containing atoms in spatially inhomogeneous external electric fields. In this context 

CN [ the subtle interplay of P-violation and time reversal (T) invariance plays an important role. If 

the entire field configuration is chosen to exhibit chirality the energies are in general shifted 
by pseudoscalar contributions which change sign under a planar reflection of the field. In 
part I we consider sudden variations of the fields and calculate P-violating energy shifts using 
Q H ' perturbation theory. In part II the adiabatic case will be treated and the connection to 

geometrical (Berry-) phases will be elucidated. To calculate the effects we use the standard 
model of elementary particle physics where the P-odd interaction arises through the exchange 
of Z-bosons between the quarks in the nucleus and the atomic electrons. We consider in detail 
hydrogen-like systems in unstable levels of principal quantum number n = 2. We study atoms 
^ with vanishing nuclear spin like ^H e+ an d with nuclear spin 1 = 1/2 like \rl. The nominal 

order of P-violating effects is 10 -5 . . . 10 -9 Hz which is determined by the mixing of the 2S\ /o 
and 2P\j2 states. However we point out that with certain configurations of the external 
fields, it is possible to enhance the P-violating energy shifts dramatically! Instead of energy 
shifts linear in the P-violation parameters we get then shifts proportional to the square root 
of these parameters. Numerically we find such energy shifts which only appear for unstable 
states to be of order 10 -5 ... 1 Hz. Under a reversal of the handedness of the external field 
configuration these P-violating shifts get multiplied by a phase factor i, i.e. the shifts in 
the real and imaginary part of the complex eigenenergies are exchanged. Application of our 
technique to hydrogen-like atoms with a nucleus of spin 1 = 1/2 yields P-violating energy 
shifts which are very sensitive to the nuclear spin dependent P-odd force, which receives 
a rather large contribution from the polarized strange quark density in polarized nuclei. 
Thus, a measurement of these energy shifts could provide an important tool to elucidate 
nuclear properties connected to the so called "spin crisis". We also present a method for 
treating degenerate perturbation theory which combines advantages of both, Kato's and 
Bloch's methods. 
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1 Introduction 



In this paper we will discuss various possibilities to obtain parity (P) violating energy 
shifts for atomic systems due to effects of the weak neutral current (NC) interactions. 
Energy shifts, i.e. frequencies can be measured in atomic physics with very high accu- 
racy. Thus P-violating energy shifts could open up a new window for measurements 
of the tiny NC effects in atoms. There are some "no go" theorems for such energy 
shifts which are related to time reversal (T) invariance. However, each theorem needs 
some assumptions, and in the following we will show how the "no go" theorems can be 
avoided by considering situations where these assumptions are not fulfilled. 



In our paper we will stay in the framework of the Standard Model (SM) of elemen- 
tary particle physics |T[ which led to detailed predictions of P-violating NC effects in 
atoms HI. By now such effects are well studied theoretically and have been verified 
experimentally for heavy atoms containing many electrons (for reviews cf. [|] ||, 
But P-violating effects have not been measured yet for hydrogen-like atoms, where the- 
oretical calculations can be pushed to very high accuracy M. Proposals for P- violation 
experiments with hydrogenic systems have been made e.g. in |8fl-|15||. In [13], a variety 



of P-violating polarization rotations was studied theoretically. It was also shown that, 
in principle, P-violating contributions to the energy levels of (unstable) hydrogen-like 
atoms in suitable external electric fields should exist. In the present work we will 
address that question in detail. 

In the days of high precision tests of the SM at high energies, e.g. at LEP and SLC 
(cf. for a review), we should list some motivations for our study: 

(1) In order to check the "running" of the SM couplings as it is predicted by 
the renormalization theory, atomic systems should be useful, providing low energy 
information complementary to high energy physics measurements. In this way a very 
sensitive test of higher order effects in the SM and a probing for new physics beyond 
the SM is possible. 

(2) P-violating effects in atoms can give interesting information on properties of 
atoms and nuclei. The experimental study of anapole moments of nuclei has recently 
become possible in this way [17], ||. P- violation in atoms may also shed light on the 
so-called "spin crisis" for the nucleons |18]. [19 . 



(3) In part II we will show how P-violating energy shifts in atoms can emerge 
as a topological effect in quantum mechanics, i.e. as P-violating geometrical (Berry-) 
phases. To study such effects should be interesting in itself. 

(4) There are many recent advances in experimental techniques of atomic physics, 
like the production and storage of "cold" heavy ion beams in storage rings (cf. |20[| 
and references therein) and the observation of interference effects with atomic beams 
(cf. ||21] for a review). Recent progress has been made on techniques of atomic mirrors 



|22|j , microtraps for ultracold atoms [f23[ and of semiconductor nano-structures, which 
could be used to realize experimental configurations as they will be discribed below. 
One could also imagine applications of the recently observed carbon nanotubes [23] or 
of the football-shaped fullerenes ||25|| . In a natural way they represent cavities with a 
diameter of a few A, which are required for the purpose of measuring our P-violation 
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effects for atoms being stored inside them. 



This paper continues the work done in [O, |T5[ and much of the formalism will be 



taken over from there. We are studying excited hydrogen-like atomic systems within 
external electric and magnetic field configurations. In the following we will briefly 



recall the methods including the Wigner-Weisskopf-approximation |26| for unstable 
states. We then examine P- violating energy shifts of atoms in external electric fields for 
two basically different experimental setups with either sudden or adiabatically smooth 
local variations of the electric field vector. Part I will be restricted to the sudden 
approximation case, while the adiabatic situation will be studied in part II which is 
in preparation. In Section ^| we consider the case of an atom inside a box with two 
segments only and show that there are no P- violating energy shifts due to T-invariance. 
In Section |3] we present results from perturbation theory for boxes with three segments. 
Our conclusions are drawn in Section [|, where we also discuss the relation of the P- 
violating energy shifts presented in our paper to those in chiral molecules. 

A short account of some of the results presented in detail here has been published 
1271. 



m 



From the methodological point of view we deal in Section |3| with degenerate per- 
turbation theory, which has been treated in the classical papers by Kato and Bloch 
|p8| , |29| and which is masterly explained in f3(|. We introduce a variant of degenerate 
perturbation theory which we find very convenient four our problems and which should 
be interesting for many other applications. Our ansatz is analogous to the introduction 
of normal coordinates in the problem of small oscillations in mechanics (cf. e.g. 



In part II, which is in preparation, we will describe P-violating effects in smoothly 
varying fields using a WKB-approximation. We relate the energy shifts of stationary 
systems via boundary conditions back to geometrical phases which the atomic wave 
functions acquire in the external field, whose topological properties in a suitable space 
are crucial. 

All our formulae are written in natural units, h = c = 1, if other units are not 
explicitly indicated. 



1.1 P-violation in hydrogen- like atomic systems in external 
fields 

We consider a hydrogen-like system, i.e. a nucleus with a single electron or muon around 
it. It is crucial for our purpose to consider the coupling of the internal motion of the 
atom, in general in an unstable state, with the cm. motion of the atom as a whole. 
Thus we will first discuss the formalism which we will use for describing this situation. 



For simplicity we will assume throughout the paper, except in Sect. |1.3| , that for both, 
the cm. motion and for the internal atomic motion, we can make the non-relativistic 
approximation. 

Our system consists of two charged particles, the lepton £ and the nucleus A with Z 
protons and TV neutrons, where I and A are bound together and experience an external 
electromagnetic field described by a 4- vector potential (A°(x), A(x)). As basic dynamic 



2 



variables we have the cm. coordinate R and the relative coordinate x: 

m e x e + m A x A 



R 



rri£ + m A 



x £ -x A , (1.1) 



where x^, x 4 are the position vectors of £, A, respectively. A simple exercise shows that 
in the non-relativistic caseQ the Hamiltonian of the system can be written as follows: 

H = ^-(P-qA(R,t)) 2 + qA°(R,t) 
2m 

+ H INT - D ■ £(R, t) - /i • B(R, t) 

+ Hq + H ra d (1.2) 

Here m and q are the total mass and charge of the atom: 

m = m A + rag 

q = (Z-l)e. (1.3) 

We denote by e (e > 0) the proton charge. The components of the momentum operator 
for the cm. motion are 

The Hamiltonian for the internal motion in the cm. system in the absence of external 
fields is denoted by H INT . The electric and magnetic dipole moment operators are D 
and /x, respectively, and £(R, t), B(R, t) are the electric and magnetic fields at the 
cm. position: 

£(R,t) = -VA°(R,t)-^A(R,t), 

B(R, t) = VxA(R,t). (1.5) 

Finally Hq 1 ^ is the Hamiltonian of the free photons and H ra d describes the coupling of 
the atom to the photons. 

In Ql. 2D we have neglected the variation of the electric and magnetic fields over the 



atomic dimensions. We have also neglected terms proportional to R x £ and R x B 
coming from the fact that the electric and magnetic fields felt by the atom are those 
in its actual rest system. Such terms give rise to the motional Stark and Zeeman 
effects. Thus for our purpose the motion of the atom should be very slow. Since we 
are interested here in the principle of P-violating energy shifts, we want to study the 
simplest possible cases. We only have to make sure that our approximations do not 
introduce any spurious sources of P-violation. 

Let us now have a closer look at the Hamiltonian of the internal atomic motion 
which we write as 

H INT = H* NT + F PV , (1.6) 

1 A detailed discussion of the cm. and relative coordinates for relativistic systems including spin 
can be found in 1331 ■ 
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where Hq nt is the P-conserving piece and H PV the P-violating one. 

In the framework of the SM the relevant effective P-violating Hamiltonian H PV is 
- to leading order — due to Z exchange between the atomic lepton and the quarks 
in the nucleus. In the non-relativistic limit one finds (cf. e.g. M; we use the notation 
as in Section 2 of |15| ): 



H PV = H$ + H<$, (1.7) 
m = ^^Q$(Z,N).{5 3 (x)(<T.p) + (<T.p)5 3 (x.)}, (1.8) 

H$ = ^=i-g^(Z,iV).{5 3 (x)(I. CT )( CT -p) + ( CT .p)(I- CT )5 3 (x)}, (1.9) 

where depends on the nucleus' spin I and H$ does not. G is Fermi's constant, 
furthermore, p and a are the momentum and spin operators of the lepton I and Q { ^ 2) 
are the weak charges of the nucleus 

Q%(Z,N) = Z(l-Asm 2 e w )-N, (1.10) 
Q%{Z,N) = -j(l-Asm 2 9 w )[Au(Z,N)-Ad(Z,N)-As(Z,N)], (1.11) 

where 9 W is the weak angle and the Aq are defined as the "total polarizations" carried 
by the quark species q = u,d, s. It is assumed in ( |1.11| ) that possible contributions of 



the c, b and t quarks in the nucleus can be neglected. We refer the reader to |l5j and 
to the references therein for more details. 

The energy levels of the hydrogen-like atom in the Coulomb approximation are 
15, 25*, 2P etc., where 2S and 2P are degenerate. These levels are split by the fine 
structure, the Lamb shift and further QED effects as well as by the hyperfine structure 
if the nuclear spin 7^0. The resulting near degeneracy of the opposite parity states 25* 
and 2P makes the n = 2 levels in principle ideal candidates for looking for P-violating 
effects. Also, the 2S state is long lived. Thus we will concentrate on the atoms in 
n = 2 states in the following. 

The relevant dimensionless parameters measuring the "wrong" parity admixture in 
the n = 2 states are 5i(Z, N), i = 1, 2, defined in (3.17) of fT5|| , which are proportional 
to Q$(Z,N). For electronic atoms we have approximately (cf. (3.21) of [p"5|j): 

6i(Z, N) ~ -Q$(Z, N) ■ 6.14 x 1(T 12 . (1.12) 

Let us discuss next the Hilbert space of state vectors for our system. We will 
consider atoms in n = 2 states, denoted \2,a), a = 1,2, ... (cf. Sect. 3.1 of Hl5|), and 



at arbitrary cm. position R, and their decay products, i.e. atoms in n = 1 states plus 
photons. A basis in our Hilbert space is then provided by the vectors 

|R; 2,a), 

|R; 1,6; t, A;), (1.13) 

where we label by k = 1, 2, ... a convenient basis set of photon states and by b = 1, 2, 
the internal n — 1 atomic states. We assume these states to be normalized to 

( R'; 2, a' | R; 2, a) = 5 (3) (R' -R)5 a , a 
( R'; 1, 6'; 7, k' \ R; 1, b;j,k) = <5 (3) (R' - R) 5 b i b 5 k >k, (1.14) 
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with all other scalar products vanishing. 

Consider now a general state vector 1 1 ) at some time t. The part of the state vector 
describing the undecayed n = 2 atoms is (cf. (3.6), (3.7) of p 



*(R,t)=X>a(R,t)|2, 



;i.!5) 



where 



V> a (R,i) = (R; 2,a|t). 



1.16) 



The wave function \1/ can be thought of as a Schrodinger wave function for the cm. 
motion, but with a number of components equal to the number N of linearly indepen- 
dent n = 2 states labeled by a. Thus we can also write \1/ as an iV-component wave 
function: 

/ Vi(R,t) \ 
$j(R, t) 

tf(R,t) < — ► . . (1.17) 

V ^7v(Rt) / 

1 states plus photons are described by multi-component wave func- 



Similarly, the n 
tions (cf. (3.8) of @) 



$ fe (R,t) = ((R; 1,6; 7, 



'1.181 



The general methods [^| for obtaining the non-hermitian mass matrix can now be 
applied (cf. Sect. 3.2 of |0|). Starting from the Hamiltonian (|1.2j ) this gives us an 
effective Schrodinger equation for the state vector \1/(R, t) as follows: 



;i.l9) 



where H e tt is a non-hermitian operator: 



H eff 

M(R, t) 
M 
D 



-(p 

2m v 



gA(R, t)f + q A°(R, t) + M(R, t), 
M - D ■ £{K, t)- n- B(R, t), 



[1.20) 
1.21) 



2,a|if /iVT 



2,a|D[2,6)), 
((2,a|/i|2,6)). 



2,6 



1.22) 



As in fL5| , quantities with a bar underneath denote matrices in the space of n = 2 
states. Thus D and ii are the corresponding matrices of the electric and magnetic 
dipole operators, respectively. In ( |1.22| ) the matrix V is the decay rate matrix for 
the n = 2 states as defined in (3.13) of fL5| . Here we have assumed that the decay 
rates from specific states | 2, a ) do not depend on their cm. motional state. This is 
an approximation which neglects e.g. relativistic time dilatation effects. But again, no 
spurious P-violation is introduced in this way. 
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In the following we will be concerned with finding solutions to ( 1.19p for various 



field configurations. In particular, we will be interested in electric and magnetic field 
configurations which are constant in time and consider the eigenstates ^(R) of H e ff 
for this situation: 

H eff ^(R) = EV(R). (1.23) 

In general the eigenvalue E is, of course, complex. The wave functions ^(R) satisfying 
( |1.23| ) correspond to states with single exponential decay. We note that a non-hermitian 



operator O does not need to be diagonalizable, i.e. in general there will be no linear op- 
erator V satisfying V~ l OV = diagonal operator. A finite iV-dimensional non-hermitian 
matrix is always diagonalizable if it has iV different eigenvalues. In the case of degen- 
eracies of eigenvalues diagonalization is in general not possible even for finite N, and 
the right (left) eigenstates form no longer a basis (cf. e.g. ||33|| ). In the following we 
will in the general analysis always assume that all non-hermitian matrices/operators 
occurring are diagonalizable. In the specific examples discussed we check explicitly 
that this assumption is valid. 

The basic question which we study in this paper is as follows: under what cir- 
cumstances can the complex energy eigenvalue E have linear terms in the P-violation 
parameters 8i j2 (Z, N). Since these parameters are very small, we will neglect terms 
quadratic in 8\^ in our calculations below. How could experimentalists look for terms 
linear in 8\^ The typical strategy which we have in mind is as follows. Consider an 
electric field arrangement £(x) and the parity-transformed one: 

£ P (x) = -£(-x). (1.24) 

Let E{£, Si, 8 2 ) be the energy eigenvalue of some level in the field £(x) and E(£ P , 81, 62) 
the eigenvalue of the corresponding level in the field £ P (x). The basic observable to 
study experimentally is the energy difference: 

AE = E{£, 8 U S 2 ) - E(£ P , Si, S 2 ). (1.25) 

Clearly, in order for AE to have a chance to be nonzero, the field configurations £(x) 
and £ P (x) must be essentially different, i.e. they should not be transformable into each 
other by a proper rotation. In other words: the field configuration £(x) must have 
handedness or chirality. 



Now we note that from considering a P-transformation we get: 

E(£,8i,8 2 ) = E{£r,-8i,-8 2 ). (1.26) 

This leads to 

£(£,0,0) = £(£ P ,0,0). (1.27) 

If we first exclude the case where there is a degeneracy of levels for <5 12 = 0, the 
eigenvalues E are differentiable in £ 1)2 at this point (cf. appendix C of [|II|) and we get: 

+ 0(8l 8i8 2 ,8l). (1.28) 

Sl,2=0 



AE = 2 ^ 



=1.2 



dE(£,6 1 ,6 2 ) 
08, 
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Thus, the measurable asymmetry ( 1.25|) is directly related to the terms in the energy 
eigenvalues E(£, 61,62) which are linear in 6\^- In practice degeneracies of levels for 
<5i,2 = do occur and have to be dealt with carefully. As we will show, such degeneracies 
can be used to obtain P-violating energy differences Ql.25| ) which are proportional to 
\f~6~i. Thus we get huge enhancements of P-violating energy shifts of order 1/ 'y/6i ~ 10 6 
for special external field arrangements. Varying the fields around these special values 
one should observe sharp resonance type phenomena. 

Of course, instead of the asymmetry (|1.25|) we could consider the asymmetry of 
the energies in the field £(x) and a field £ R (x) obtained from £ by reflection on some 
plane. Maybe it is also possible to make a continuous change from £ (x) to £ P (x) or 
£ R (x) and observe the corresponding continuous change of the energy eigenvalues. 



1.2 Atoms in a constant electric field 

Let us consider a simple case to start with: an atom in a constant electric field £. We 
assume the atom to be confined inside a region B which can be of arbitrary shape. 
This confinement of the atom in B could, for instance, be realized by having the walls 
of a suitable trap containing the atom built as "mirrors" (cf. ||22|| ). We represent the 
"vessel" or "trap" containing the atom by a potential V(,(R) which is zero inside the 
vessel and infinite outside. This leads to the following boundary condition for the wave 
function ^(R): 

(R) = for R G dB (1.29) 

where dB denotes the boundary of B. In the following ( |1.29| ) is always understood to 
be imposed. The effective Hamiltonian ( |1.20| ) reads now 

p2 

H eff = — ~qR£ + M -B-£. (1.30 

2m 

The eigenvalue equation (|1.23|) can easily be discussed for this case. Clearly, we can 
diagonalize the operators governing the cm. motion, 

ffE " = ^-« R ' £ - < L31) 

and the internal motion, 

M = Mo -!!■£, (1.32) 

separately, form the tensor product of the corresponding eigenvectors and add the corre- 
sponding eigenvalues. Indeed, let n (R), n = l,2,...,bea complete set of orthonormal 
eigenvectors of H EXT : 

H EXT un) = K XT MK)i (1-33) 

where E EXT are the eigenvalues which are, of course, real since H EXT is a hermitian 
operator. Let \(3) be the right eigenvectors of AA to eigenvalues E INT (f3) 1 which can be 
complex: 

M = E INT (P) |/5), (/? = 1, 2, N). (1.34) 
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Then the eigenvectors \I/ of H e ff ( |1.30[ ), ( |1.23| ) are 

* /9in (R) = ft (R)|/9) = |^n), (1.35) 

with eigenvalues 

i?(/?,n) = C T + ^ T ^ 1 (1-36) 

i.e. we have 

tf ef * Aft (R) = E(P, n) ^, n (R). (1.37) 

Here we can assume without loss of generality that the electric field points in 3- 
direction: 

£ = £e 3 . (1.38) 
The eigenvalues and eigenvectors of the corresponding mass matrix 

M=M -D 3 £ (1.39) 



were discussed extensively in JT5[. For later use and to fix the notations we recall the 
main features here. 

For an electric field ( |1.38j ) the mass matrix A4 ( |1.39| ) commutes with the 3- 
component F 3 of the total angular momentum operator: 

[M,F 3 ]=0. (1.40) 

We can diagonalize A4 and F 3 simultaneously. The right eigenvectors \oc, F 3 ,£e 3 ) 
satisfy: 

F 3 \a,F 3 ,Ee 3 ) = F 3 \a, F 3 ,£e 3 ), 

M\a,F 3 ,£e 3 ) = E INT (a,F 3 ,£e 3 )\a,F 3 ,£e 3 ); (1.41) 
the left eigenvectors (a, F 3 ,£e 3 \ satisfy: 

{a,F^£e 3 \F 3 = (a,F^£e 3 \F 3 , 

(a,F^£e 3 \M = {a ) F 3 ^£e 3 \E INT {a,F 3) £e 3 ). (1.42) 

Here we label the eigenvectors by (a,F 3 ) instead of j3 as above, with a = 1,2,... 
numbering the eigenvectors to the same F 3 . We will assume that there is no degeneracy 
of the complex energy levels for fixed F 3 : 

E INT (a, F 3 , £e 3 ) ^ E INT (a', F 3 , £e 3 ) for a ^ a'. (1.43) 

This is true for all systems in external fields which we will consider. The condition 
( |1.43j ) guarantees that the right eigenvectors ( |1.41| ) (left eigenvectors ( |1.42|) ) form a 
basis (the dual basis) in the n = 2 subspace where the normalization can be chosen 
such that 

(a, F^£e 3 \a', F 3 , £e 3 ) = 5 a ^ 5 F3>F >. (1.44) 
We introduce quasiprojectors (cf. (3.28) of 



P(a,F 3 ,£e 3 ) = \a, F 3 ,£e 3 )(a, F 3 ,£e 3 \ (1.45) 
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which satisfy 



E(a,F 3 ,£e 3 )E(a',F^£e 3 ) = 8 a , a , 8 F ^ n W{a, F 3 , £e 3 ), 

Tr P(a,F 3 ,£e 3 ) = 1, (1.46) 

but are not necessarily hermitian. We have 

1 = ^P(a,F 3 ,£e 3 ), (1.47) 

a,F 3 

M = F 3 ,£e 3 )P>,F 3 ,£e 3 ), (1.48) 

a,F 3 

E 3 = £F 3 P(a,F 3 ,£e 3 ). (1.49) 

a,F 3 

The total energies (|1.36|) are 



E(a,F 3 ,£e 3 ,n) = E° XT + E INT (a,F 3 ,£e 3 ). (1.50) 



It was shown explicitly in Sect. 3.4 of |15] that the complex energy eigenvalues 
E INT (a, F 3 ,£e 3 ) have no contributions linear in the P-violation parameters 5i^- This 
was derived by considering a reflection R on the 1-3-plane for the atom at rest. Now 
we have the atom in a box which need not be R- invariant. But for the constant 
electric field we have no coupling between the cm. and the internal motions in our 
approximations. Thus even if the external field arrangement is not R-symmetric here, 
we still get no terms linear in <5 12 in the total energy eigenvalues E(a, F 3 , £e 3 , n) (|1.50|) 
for the constant external field. 

In the following we will investigate situations where a non-trivial coupling between 
cm. and internal motions is present. We will see that indeed for suitable external field 
arrangements having no reflection symmetry the total energy eigenvalues can get terms 
linear in the P-violation parameters 5i 2 . 



1.3 General consequences of T-invariance 

In this Section we will work out some general consequences of T-invariance for sta- 
ble and unstable atomic systems in a static external electromagnetic field, described 
by the 4-potential (y4^(x)) = (A(x), y4°(x)). Let us start by considering the atomic 
system coupled to the radiation field at the level of quantum field theory (QFT). The 
Hamiltonian describing this system is 

H[A] = H (t) + H%A], (1.51) 

where Ho, which includes the P-violating term H PV , is the Hamiltonian for = 0, 
and H' is the interaction term with the external field. In terms of the fundamental 
lepton and quark fields the latter reads 

H%A] = J d 3 xj M (x,t)^(x), (1.52) 
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with jfj,(x, t) the electromagnetic current operator. For static fields the total Hamilto- 
nian is time-independent, as we know from Noether's theorems. In the following we 
will work in the Schodinger picture of quantum mechanics, identifying the Hamiltoni- 
ans there with those of (|L51| ) at t = 0. Thus we set H = H (0) and H'[A] = H'[0, A}. 
We will now assume that another splitting of H can be made: 

H[A] = H Q [A] + H rad , (1.53) 

where £/o|y4] describes stable atomic systems in the external field A^(x.) and H ra d their 
coupling to the radiation field. This coupling will render the above systems except the 
lowest level unstable. Let | (o), a, A) be the eigenstates of i?o[A] : 

H [A) | ip), a, A) = E (a,A) \ (o),a,A), (a = 1,2,...), (1.54) 

where E (a, A) are real eigenvalues. In practice it turns out to be convenient to work 
with some suitable basis of orthonormal states 

\a,h) (a = l,2,...;n = 1,2,...) (1.55) 



for the eigenstates in ( |1.54j ). Below, the index a will be chosen to be related to the 



internal state of motion of the atom, n to the state of the cm. motion. We denote by 
Ti. the Hilbert space spanned by the states ( |1.55| ). 



With the inclusion of H ra( i of (|1.53| ) the states | a, h) become in general unstable. 
Using the standard Wigner-Weisskopf procedure (cf. Sect. 3.2 of |15[]) we find the 
effective Hamiltonian H e g describing the time evolution of the undecayed states \a, n ) 
through the (approximate) solution of Schrodinger's equation with the initial condition 



of having some combination of the states | a, h) at time t = (cf. (3.9) of [[15 



i^\t) = H[A]\t), 

l t = °) = ^Zc a) n\a,h). (1.56) 

u.h 

For t > but neither too small nor too large times we get 

( a, h | exp(— iH[A] ■ t) \ a', n' ) ~ ( a, n \ exp(— iH e ff[A\ ■ t) \ a', h'). (1-57) 

This represents the general definition of the effective Hamiltonian, which reduces to 
( |1.20| ) in the non-relativistic limit. 

We consider now time reversal (T) invariance. From QFT we learn that there exists 
an antiunitary operator V(T) which transforms the Hamiltonians as follows: 



V{T)H'[A] V- 1 (T)) t = H'[A T ], (1.58) 
(v(T)iJ[A]V- 1 (T)) t = H[A T ], (1.59) 
V(T)^ [^]^ 1 (T)) t = H [A T }. (1.60) 
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Here A T denotes the T-transformed potential: 

A°(x) = l°(x) 



A T fx) 



-Afx). 



:i.er 



From ( 1.54]) and ( |1.60| ) we find that V(T) transforms the eigenstates | (o), a, A) of H [A] 
into eigenstates | (o), a, A T ) of ff [^4 T ]- Recall that E (a, A) are the eigenvalues of the 
hermitian operator i?o[^4] and thus real. We assume now that the space 7i spanned by 
the basis vectors \a,n) in ( |1.55|) is large enough to contain both the states | (o), a, A ) 
and | (o), a, A T ). The action of V(T) o:n the states | a, n ) with chosen phase convention 
is then given by some matrix T over 7i: 



V(T) \a,n) = | a' : n'}f a 'h',o 
The antiunitarity of V(T) implies: 

f f T = 1. 



;i.62) 



1.63) 



The effective Hamiltonians H e ff[A] ( |1.57| ) and H e jf [A T ] can be identified with non- 
hermitian matrices over Ti 

(1.64) 



H eff [A] = [{a,n\H eff [A]\a',h') 

From ( |1.57| ), (|1.59| ) and ( |1.61|) we now get as a consequence of T-invariance: 

f^H eff [A T ]f) T = H eff [A}. (1.65) 

This is the generalization of (3.42) of [O] for atomic systems in arbitrary static external 
fields. 

Let now {E(a,A)} be the spectrum of eigenvalues of H e jf[A] and P(a,A) the 
quasiprojectors corresponding to E(a, A), a = 1, 2, .... 



He ff [A] = J2E(a,A)W(a,A). 



;i.66) 



An argument using the resolvents of H e ff[A] and H e ff[A T ] analogous to (3.43) ff. of flBH 
shows that with a suitable numbering we have 



E(a, A^ 
(f^(a,A T )f 



E{a,A), 
W(a,A). 



;i.67) 
:i.68) 



These equations embody the general consequences of T-invariance for the complex 
(real) eigenvalues and quasiprojectors (projectors) of unstable (stable) atomic states 
in static external fields. 

Let us consider next a general improper rotation R around some point x (0) , which 
we choose without loss of generality as x (0) = 0: 



R : x 

det R 



-1. 



1.69) 
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For the external field we set: 



<(x) = A^sr 1 *), 

A R (x) = i2 • A(i2 _1 • x) (1.70) 



The QFT-operator U(R) corresponding to R transforms the Hamiltonian (|1.51 ) as 
follows: 

U(R) H[A, 81, 8 2 ] U~ 1 (R) = H[A n , —81, —82], (1.71) 

where we indicate now explicitly the dependence on the P- violation parameters 8^2- 
Assuming the space 7i to be large enough to contain for a state also its U (R)-transform 
we can represent U(R) in 7i by a matrix 1Z\ 

U(R)\a,h) = \a',n')TZ a ffi>,ah, 

K ] K = 1. (1.72) 

For the effective Hamiltonian this leads to 

T$ H eff [A R , -81, -82} K = H eff [A, 6 U 8 2 }. (1.73) 

Considering again the resolvents of H e ff[A, Si, 8 2 ] and H e g [A R , — Si, — 8 2 ] we find that 
the numbering of the spectra of complex energy eigenvalues can be chosen such that 

E(a,A n , -81,-821 = E(a,A,8i,8 2 ), (1.74) 
-R)V{a,A n , -81, -82)K = W(a,A,5 u 8 2 ). (1.75) 



From ( |1.67 ) and ( 1 . 74 ) we find the following 



Theorem: If the external field A has the property that the T-transformation for it is 
equivalent to some improper rotation R 

A T = A n , (1.76) 

then 

E(a, A, 81, 82) = E{a, A, -8 X , -8 2 ). (1.77) 

For the case that the energy level E(a, A, 81, 8 2 ) is not degenerate for 8\ = 8 2 = 0, this 
implies vanishing linear terms in 8^2' 



d 

—E(a,A,8 u 8 2 ] 



= 0, (j = l,2). (1.78) 

<5i=<5 2 =0 



As an example for the use of this theorem we consider a hydrogen atom at rest at 
x = in static and spatially homogenous electric and magnetic fields S, B. Let R be 
the reflection on the plane through x = spanned by £ and B. We can write 

A°(x) = -x£, 

A(x) = -ixxB, (1.79) 
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which fulfils 

^(x) = A£(x). (1.80) 

Furthermore, for £ ^ and B ^ there are in general no degeneracies of levels for 
Si — S 2 — 0. Thus we conclude that the energy eigenvalues of the stable and unstable 
levels of the hydrogen atom in a constant electric plus magnetic field have no terms 
linear in the P- violation parameters . 

To illustrate the general T- and R-transformations we will now close this Section 
by discussing these for the non-relativistic effective Hamiltonian ( |1.20|) . The space 7i 



(cf. (|1.55|) ff.) is then spanned by the vectors |R;2, a) of ( |1.13|) . We can choose a 



convenient orthonormal basis set of real functions £n(R) for the cm. motion and set: 

|a,n)=&(R)|2,a). (1.81) 
The T-transformation Ql.621) is represented by 



V(T)\a,h)=&(R)\2,a')T a , a , (1.82) 

where T is the matrix describing the T-transformation of the n = 2 atomic states (cf. 
(3.33) of fl5fl). Thus we find here 

Xi'h',an = Sn'nUa'a- (1.83) 

It is now easy to verify ( |1.65| ) explicitly for H e ff[A] of ( |1.20| ). Thus, as it should be, the 
general results (|1.67|) , (|1.68|) hold in the non-relativistic limit. 

Let us furthermore define the matrix TZ of the R-transformation ( 1.69| ) for the n = 2 

states 

U(R)\2,a) = \2,a')E a , a . (1.84) 

and (7%,-) by 

^(i?" 1 -R)=^(R)^k- (1-85) 
Then the matrix 1Z of ( |1.72 ) is given by 

Tla'h',ah = Tl'n'h'Ra'a- (1.86) 



This is the matrix which is to be used in (|1.74|) and ( |1 . 75|) in the non-relativistic limit. 



In the following we will frequently make use of the general consequences of the T- 
and R-transformations spelled out in this Section. 



2 Energy shifts for piecewise constant electric 
fields; exact calculation 

We will now consider simple electric field arrangements which are neither parity- nor 
reflection-symmetric. We imagine a rectangular box which we subdivide into several 
segments, in each of which the £ -vector is pointing in a fixed direction, with its modu- 
lus also having in general different values. Within any of the segments the electric field 
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is assumed to be homogeneous, which results in a sudden variation at the interface, 
where the segments are connected to each other, and at the box surface which forms 
the boundary of the system. Of course, Maxwell's equations require then charge dis- 
tributions on the interfaces. Clearly it would be difficult to realize this experimentally, 
and we choose this situation only since it is simple to be analysed from the theoretical 
point of view. A more realistic experimental situation could be taken into account 
by including corrections to the sudden variation by means of perturbation theory in 
a small parameter describing e.g. the extension of the transition regions between the 
segments in relation to the total size of the segments. 

Let the lengths of the edges of our rectangular box B be ^1,2,3 in the 1,2,3-direction 
(Fig. 1): 

box B : < Rj < Aj, (j = 1, 2, 3). (2.1) 
In 1-direction we subdivide the box into segments numbered by a: 

segment S M : R ( {~ 1] < R x < R[ a) , < Rj < Aj, (j = 2, 3), (2.2) 

where a = 1, ...,cr max and 

R m = < Rf } < ... < R[° max) = A x . (2.3) 

This is illustrated in Fig. 1 for a max = 3. The electric field in S <a) is denoted by £ {a) . 

We will consider both neutral atoms and charged ions in such segmental boxes. The 
effective Hamiltonian ( |1.30|) reads now 

H eff = -^P 2 + M - (qR + D) • £ 6(^ CT) - Rj e(R 1 - R^), (2.4) 
2m a 

The problem is to solve the eigenvalue equation ( |1.23| ) with the boundary condition 
( |1.29| ). There is a direct way to solve this problem by making an ansatz of linear 
combinations of constant field wave functions in each segment. Continuity conditions 
to be imposed on the wave function and its first derivative at the interfaces lead to a 
determinant equation from which the eigenenergies may be extracted. This technique 
will be outlined in the next paragraph and yields the minimum number of different 
segments, i.e. electric fields required to obtain P-odd effects. We will show that for 
boxes with 2 segments no energy shifts linear in the 5i (i = 1, 2) occur in our theoretical 
framework. However, for boxes with 3 segments such energy shifts do occur as we will 
prove in Sect. [3] using perturbation theory. 

2.1 A neutral atom in a box with 2 segments 

In this Section we consider a neutral atom in internal states of principal quantum 
number n = 2 in a box with two segments. The effective Hamiltonian ( |2.4j ) reads now: 

H eff = ^P 2 + M - D • S(R), (2.5) 

where 

£(R) = £ (1, e« - Ri) e{R t ) + e (2) e(A 1 - e{Rt - r[ 1} ). (2.6) 
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We have to solve ( OB ) with the boundary condition ( 1.29| ): 

H eff V(R) 



VCR) 



EV(R), 
for R e dB. 



(2.7) 
(2.8) 



The Hamiltonian ( p.5| ) commutes with the part describing the motion in the 2- and 
3-directions: 



(2.9) 



Here and in the following the boundary condition ( |2.8|) is always understood to be 
imposed. Thus we can diagonalize H e ff and H 23 simultaneously. The normalized ei- 
genfunctions of H 23 are 



( Pn 2 ,n 3 (R2, R3 

(n 2 = 1,2,3, n 3 = 1,2,3, ...). 



We have 



. n 2 TcR 2 . n 3 7fR 3 
sin — ; • sin 



An 



A, 



(2.10) 



#2,3 (R 2 , R 3 ) 

E 2t3 (n 2 ,n 3 ) 



E 2 , 3 (n 2 ,n 3 ) 4> n2in3 (R 2 , R, 
fn 3 ir\ 2 
\A 3 J 



3), 



1 

2m 



n 2 n 
~A7 



+ 



(2.11) 



We can now make a product ansatz for the eigenf unctions of H e ff 
*(R) = 

<P(Ri) = 



(R^R^^Ri), 
f (1) (^i) for < Ri < Rf\ 



U2) 



(Ri) for R[ 1} <Ri<A h 



f2.12) 



where <fi(Ri) has N components (cf. ( |1.17| )), corresponding to the number of linearly 
independent n = 2 internal states. 

Let us denote by |0),/?) the right eigenvectors of the mass matrices in the Sec- 
tors cr = 1,2; corresponding to the internal, in general complex energy eigenvalues 

E INT ((v), (3): 

(M - D • £ (<r) ) |K/?) = E INT ((a),(3) \(a),(3), ((3 = 1, ...,N). (2.13) 

The corresponding left eigenvectors are denoted by (0),7|, (7 = 1, N). We have 

(<Zyrf\(r),P) = 6 7 p for a = 1,2 (2.14) 

(no summation over a). Furthermore we define matrices U (1 ' 2) and U {2,1) with matrix 
elements: 



r(2,i) 



U^(5 U S 2 ) 



(«,7|( 2 ),/3), 



(2.15) 



We indicate explicitly the dependence of these matrices on the P-violation parameters 
5i 5 2- Of course, for the case of nuclear spin 1 = there is no dependence on 5 2 . As we 
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discussed in Sect. [L72|, the complex energies E INT ((&), j3) have no terms linear in 8± j2 . 
We make now the following ansatz for (1,2) (i?i) in ( p. 12 ): 



> (1 W 



EC, 



a) 



sm 



'P 

1(2) 



SCfsin Qf{R x -A l ) |( 2)) (3), 



(2.16) 



where Qp and ^ are constants to be determined. With ( |2.12| ) and ( |2.16| ) we already 
satisfy the boundary condition ( |2.8| ). Let us now insert the ansatz (|2.12| ), ( |2.16| ) in 



and let E be the sought eigenvalue of H e g. We find that the Qg in ( p,16| ) are to 
be determined from 



E-E 2i3 (n 2 ,n 3 )-E"' T (W),I3), (a = 1,2; /3 = 1, N). (2.17) 

For definiteness, let Qf be the root of (|2.17| ) with positive real part and if this is 
zero, the one with positive imaginary part. Finally, the constants Cg in ( |2.16| ) are 
determined by the transition conditions at the interface R\ = R± . There ^(R) and 
d^/(R)/dR\ = di^(Tl) must be continuous. These conditions lead to: 

£ sin [QfR^] |(D, P) = ECsin^^-^)] |(2),/5), (2.18) 



P 



QfR?]\(D,P) = CpQf cos [Qf ( R i ~ A i)] K 2 ), P)- (2-19) 



P P 

Multiplying ( |2.18| ) and ( |2.19| ) with ((2), 7] and ((1), 7I, respectively, where 7 = 1 
we get the following matrix equation: 

M e (E,5 1 ,5 2 )-C = 0, 

where 



,...,N, 



(2.20) 



M C {E,6 1 ,S 2 ) 



U^(5 1 ,S 2 )-S m (E,5 1 ,S 2 ) -S™(E,5 U 5 2 ) 

C W (E,5 U 5 2 ) -U^{6 1 ,6 2 )-CV{E,5 1 ,6 2 ) 



S — \ S,l 5 •••) SAT 5 S>1 1 •••S>. 



(1) A(2) 



-(2) 

N 



(2.21) 
(2.22) 



In ( p.21| ) and C* (<7) (a = 1,2) are diagonal N x N matrices, with the diagonal 
elements given by 



S ( ^(E,5 1 ,5 2 ) 
C^{E,5 X ,5 2 ) 



sm 



sm 



V/3 K l 



Qf - A, 



Qf cos QfR? 



C^(E,8 1 ,8 2 ) = Qf cos [Qf (r? - A, 



(2.23) 
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Note that in the matrices S i<y) , C {a) the dependence on the parameters 8^2 enters 



only through the internal energy eigenvalues E INT ({v), (3) (cf. ( 2.17)) ), which are even 
under the replacement 8^2 — > — ^1,2- Therefore we have 



S^\E,6 1 ,6 2 ) = S<'\E,-6 u -6 2 ), 
C^(E,5 1 ,5 2 ) = C^{E,-8 u -8 2 ), (cr = 1,2). 

The condition for ( |2.20|) having a nontrivial solution-vector £ is: 

det M C (E,8 U 8 2 ) = 0. 

The roots of this equation give the complex energy eigenvalues of H e ff 
to see if these roots have terms linear in 81 2 ■ 



(2.24) 

(2.25) 
51). We want 



2.1.1 Results for neutral atoms with 1 = 

Let us first consider the case of a neutral hydrogen-like atom with vanishing nuclear 
spin. This is not realizable for conventional atoms but an atom consisting of a tt + 
meson and an electron is a physical system to which our considerations apply. For this 
case the representations of A4 and D may be found in Appendix B of [12]. They are 
8x8-matrices and thus the continuity conditions for the wave function and its derivative 
w.r.t. Ri lead to a set of 16 linear equations (cf. ( |2.20| )): 



From the determinant condition 



detM c (#A) = 



(2.26) 



(2.27) 



we get the energy eigenvalues. As we show in Appendix [TJ we can simplify the structure 
of M c considerably if we choose a suitable set of basis states having simple transforma- 
tion properties under T (time reversal) and R (reflection on the 1-3-plane). We show 
in Appendix |B| that ( 2.27Q is equivalent to an equation having no linear terms in Si. 



det M' e (E, 81) = f (E) + 8 2 J 2 (E) + 0(5$) = 



(2.28) 



Now we would like to use the implicit function theorem to calculate from ( |2.28| ) the 
dependence of the solutions E on 8x- 

However, we have to be careful. It turns out that for 81 = the energy levels are 
twofold degenerate, i.e. fo(E) has double zeroes. These degeneracies are even preserved 
to order 81. We show all this in Appendix B, where we use the T-symmetry to derive 
that det M' c = is equivalent to 



det M« sym (E, 81) = 0, 



^2.29) 



where M^ sym = —(M^ sym ) T is an antisymmetric matrix. Now, from a basic theorem 
of linear algebra (cf. Appendix 0), the determinant of an antisymmetric matrix with 
even number of rows is always a square: 



det Mc Sym (E, 81) = \g(E,8i 



[2.30) 
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Therefore the solutions E of ( |2.29| ) are twofold degenerate for any choice of £ (1) ,£ (2) . 
Furthermore we find that 



g(E,5 1 )=g (E)+5 2 1 g 2 (E) + O(8t 



(2.31) 



with go{E) having only nondegenerate zeroes. The implicit function theorem implies 
then that the zeroes of g(E,5i) have no terms linear in Si: 



dE(6 1 ) 



dSi 



0. 



(2.32) 



<5i=0 



To conclude: Using the R and T symmetry relations we find that, for a box with two 
electric fields £ (1) ,£ (2) , the complex eigenenergies of the atom have no terms linear in 



2.1.2 Results for neutral atoms with 1 = 1/2 

In this subsection we will consider the case of a neutral atom with / = 1/2, as e.g. \B.. 
The eigenenergies are again obtained from the condition 

det M C (E,5 1 ,6 2 ) = (2.33) 

for the continuity matrix M c given in fl2.21| ). Now M c is a 32x32 matrix. Using the ± 
eigenstates as defined in Appendix ^ and their T and R symmetry relations we show 
in Appendix [B] that ( |2.33| ) is equivalent to 

det M' C (E, 5 U 5 2 ) = f (E) + ^ 5 i 5 j f 2 , ij {E) + terms of 4th order in 8 h2 . (2.34) 

hi 

Here fo(E) has (in general) no double zeroes and thus we can immediately conclude 
from the implicit function theorem that the energy eigenvalues, i.e. the solutions of 
( |2.34j ), have no terms linear in 5\ y2 : 



dE{Sx,S 2 ) 



d8 j 



= 0, (2.35) 

<5i=5 2 =0 



Thus, within our theoretical framework, we find that, for a jH-atom in a box with 2 
segments only there are no energy shifts linear in 5\ 2 . 



2.2 A charged ion with nuclear spin I = in a box with two 
segments 

In the real world conventional one-electron atoms with vanishing nuclear spin are always 
electrically charged. Therefore we investigate now the modifications of the results of 
Section |2.1.1| due to the cm. motion being influenced by the external electric field. 

Let q = e(Z — 1) be the charge of the ion. The effective Hamiltonian is given in 
( \T4y where we choose 2 segments, i.e. £(R) as in ( p.6|) . In Sect, ^j] the cm. part of 
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the wave functions was built up by sines and cosines. Now, for a charged ion, we need 
the wave functions for a potential linear in the coordinate. 

Let us consider first the general one- dimensional Schrodinger equation with a po- 
tential linear in the coordinate, 



d 2 



OR 2 



-2m(qSR+ E) 



ip{R) = 0. 



^2.36) 



The solutions of ( |2.36|) are obtained with the help of the Airy functions defined by 

m(z) = yi (i_ 1/3 (c)-wc)), 



Bi(z) 



I-l/s(C)+W0 



^2.37) 



for z > and 



M(Z) = |^(j-l/3(C) + ^ + l/3(C)) 

Bi(z) = 



-I j_ 1/3 (0- j +1/3 (0 



for z < 0, where J p (I p ) denote the (modified) Bessel functions, 



2c I ^ 1 3/2 



2mq£. 



The general solution of ( |2.36| ) is then 

i>(R) = Z+MR, E, £) - {R, E, £), 

where 

ip + (R,E,£) = Ai(R + 2mE/c 2 ), 
ip„(R,E,£) = Bi(R + 2mE/c 2 ), 



(2.38) 



^2.39) 



(2.40) 



(2.41) 



and £-t are constants. If we impose as boundary condition the vanishing of the wave 
function at R = and R = A, 



V>(0) = j>{A) = 0, 
we get corresponding eigenenergies E by solving the equation 



det 



V+(o,£,£) il>-{o,E,e) 

MAE,£) il>-(A,E, £) 



0. 



(2.42) 



(2.43) 



Let us number these eigenenergies by E m (A,£) and the corresponding normalized 
eigenfunctions ( |2.40|) by <p m (R,£), where m = 1,2,.... We can choose the functions 
4> m (R, A, £) to be real and normalized: 



{m,A,E\m',A,£)= / dR <f> m (R, A, £) <j> m >(R, A, £) = S„ 
Jo 



(2.44) 
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Now we come to our problem of finding the eigenvalues of H e g ( |2.4| ) for a box with 
2 segments. We make the following ansatz for the complete wave function \1/(R): 

9(R)-l * (1)(R) f ° r °< R ^< R i^ (2 45) 

V W " i *W(R) for R[ 1} <R,< A lt {2Ab} 



^ (CT) (R) 



. , ». j ,m 3 ^mj (-^2 , ^2 , ^2 CT) ) 0m 3 (i?3 , ^3 , £{ ) 



f3,m2,m 3 



(°0 

+,/3,m 2 ,m 3 1 



^ = 1,2), 



u(o-) 



^(Rx,E^% 2>m3 ,£n\\(a),(3), 

(2.46) 



where we use the same notation as in Sect. j2.1| , and we have with £7 the sought energy 
eigenvalue of H e ff\ 



Furthermore, we set 



E mz {A^£^)-E' NT ^)^). 



(2.47) 



>±,/3,7772,m 3 



•0^(0, E^ m2tm3 , £[ \ 



This ansatz fulfils already the boundary conditions at the complete surface of the box. 
We still have to impose the continuity of ^(R) and d^f(H)/dRx at the interface, i.e. 
at Rx = Rx x) . This gives conditions analogous to ( [2. 18] ), ( |2.19| ), which we can write as 
linear equation for the now infinite component vector (: 



where 

(3,(3' 

Explicitly we have: 

M c (E,5x) 



M c (E,6x)-( = 0, 



c = /Y (1) c {2) 

^> \^/3,m 2 ,m 3 > ^>/3' ,m' 2 ,m' 3 

, N; m 2 , m' 2 = 1, 2, m 3 , m' 3 = 1, 2, 



U^iSx^iE^x) -S (2) (E,5x) 
C w {E,5x) -U ih2) {6x)G (2) {E,6x) 



(2.49) 



(2.50) 



(2.51) 



U^(5x) = ((H,a|M,/?)(m 2 ,A 2 ,^V2,^2,^ ) )(^3,A3,^ CT V / 3^3,^ T) )) , (2-52) 



S (a) (E, Sx 



>af3 



/3,777 2 ,777 3 ' ? /'+(-^l ) ' ^ftm 2 ,m 3 i ^ 



-c ) / 3, m2 , m3 V'-(^i ,e, 



/3,7772,771 3 ' 



(2.53) 







>a3 



dRx 

) 

/3,7772,777 3 



£+ ) ,P,rn2,rn 3 ' l P+(R'l-' ^f3,m 2 ,m 3 ' ^1 ') 



-e (CT) 



Ri=R m ^m2,7n^m 3 ,m^ J • (2.54) 
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Now we choose as in Sect. |2.1.1| and Appendix [B] the ± basis states ( |B.6| ) for the 



internal motion of the ion. The matrices |7 (<T ' T) can then be split into submatrices 
jj{",t)ts ^ s g |_|_^ _ j-j ag j n (p^ ^ (|B.8|) . Of course, all these matrices are now infinite 



dimensional. It is easy to check that the matrices JJ (,7,T)rs satisfy again the relations 
( |B.10|) -( |B.12 ). Thus for any finite truncation in the transverse modes, say 



m 2 ,m' 2 < K 2 , 

m 3 ,m' 3 <K 3 , (2.55) 

we can immediately take over the methods of Sect. |2.1.1| and conclude that there are 
no shifts in the eigenenergies E which are linear in 8\. 

Of course, we would now like to take the limit K 23 —>■ oo. We will not give a rigorous 
mathematical treatment of this limit but discuss it in physical terms. Suppose we have 
an eigenenergy E. Then we see from ( [2 .47 ) that the energies Ej£ m ma will be more and 



more negative for higher and higher transverse energies E m2 , E m3 , i.e. for 777.2, m 3 —> 00. 
But for large negative Eg^^ the Airy functions ( |2.37|) , (|2.38|) will be exponentially 



decreasing or increasing. Clearly, the corresponding contribution in ^ (<7) (R) ( |2.46| ) 
must in essence be exponentially decreasing when going away from the interface or 
surface of the box, and this decrease will be the faster the more negative Ep mam3 is. 
We can write the energy eigenvalue E as 

E = X f \ ' . 2.56 

Choosing now some large enough K 2 , K%, we see that the contribution from the high 
transverse modes with m 2 > K 2 , m 3 > K 3 will in essence come from the interface 
and/or the surfaces and can be made arbitrarily small in the integrals (|2.56[ ). Here it 
is important that H e g contains no parts with high frequency in any direction. Thus 
we estimate that the limit K 2 ^ — > 00 should be harmless. 

In conclusion: In this Section we found that, within the framework of our calcula- 
tions, there are no energy shifts linear in the P- violation parameter 8\ also for a charged 
ion with nuclear spin / = in a box with two segments. 



Energy shifts from perturbation theory for atoms 
in a box with three segments 



In Sections [O] and fX^ we considered boxes with two segments, but did not find any 
P-violating energy shifts. This is somewhat surprising, since the setup with two electric 
field vectors together with the fixed direction marked by the external geometry does 
form a chiral situation. But apparently, this is not sufficient to have energy shifts 
revealing parity violation. 

The above results require another segment to be added and to consider a situation 
where one can form a nontrivial trihedral out of the electric field vectors alone. For the 
simplest case of neutral atoms with 1 = 0, this adds another 16 rows and columns to 
M r ( p.21 ), (|B.13| ), as now continuity conditions have to be imposed at two interfaces 
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and the wave function 0f ) of the middle segment consists of sine and cosine terms. 
The result is a determinant equation for a nontrivial 32x32-matrix, which we have not 
analysed, since it is rather complicated. Instead, we turn in this Section to simpler 
calculations based on Rayleigh-Schrodinger perturbation theory. 

Let us consider a box with three segments (Fig. 1) where the electric fields are 

= £e 3 + £' (CT) , (a = 1,2,3). (3.1) 

We assume that the electric field in 3-direction, £e 3 is dominant: 

|£ V) | < \£\, ((7 = 1,2,3). (3.2) 

We will then write down a perturbation expansion with respect to the terms induced 
in H eff ( FH) by the fields S' {a) . Indeed, let us write H eff (pi]) as follows: 

Heff = H e ff,0 + H efftl , (3.3) 

where 

H efffi = -^P 2 + M - (qR 3 + D 3 )£, (3.4) 

# e#)1 = -(<zR+ D) • ©(^i CT) - ^i) ©(^i - Ri~ l) )- (3-5) 

<T=\ 

Denoting by R the reflection on the 1-3-plane x<i = A 2 /2, we see that H e ff^ is R 
symmetric whereas H e ff t \ in general is not. Thus we expect from the results of Sect. 
|1.3| that the eigenenergies of H e ff will not have linear terms in 5 li2 but that such terms 
will be generated upon inclusion of H e ff i with a suitable arrangement of fields £ (<T) . 

The idea is now to calculate 5j-lmear contributions to the eigenenergies E of H e ff 
in perturbation theory based on the separation ( |3.3| ). The eigenvectors as well as 
the eigenvalues E m of H e ffQ may be determined numerically in a simple way and the 
series of higher corrections to E {0> contains just matrix elements of # e ff,i w.r.t. these 
eigenvectors. 



3.1 Results for a neutral atom with 1 = 

Let us again consider first a neutral atom with nuclear spin 1 = 0. The eigenvectors 



of Hefffi are then easily written down, using the results of Sect. |L2 

\a,F 3 ,£e 3 ,n) : = 



J n (R)|a,F 3 ,£e 3 ), 



fR) 



i=i 



The corresponding eigenvalues E (0) (a, F 3 ,£e 3 ,n = (ni,n 2 ,n 3 )) satisfy 



H e ff,o\a,F 3 ,£e 3 ,n) 
E^(a,F 3 ,£e 3 ,n) 



F (0) (a,F 3 ,£e 3 ,n)|a,F 3 ,£e 3 ,n) 



1 JL /n,-7r x 2 



2m 



E 



A, 



+ E 1NT (a,F 3 ,£e 3 ). 



(3.6) 

(3.7) 
(3.8) 
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We know from the results of Sections 3.3 and 3.4 of \TE\ that E INT (a, F 3 , £e 3 ) has 
no terms linear in 5\ and that these energy eigenvalues for F 3 and — F 3 are degenerate. 
These statements are then also true for E (0) (a, F 3 , £e 3 , n), i.e. we have 

E m (a,F 3 ,£e 3 ,n) = E m (a, -F 3 ,£e 3 ,n), (3.9) 

and, neglecting terms of order Sf, E (0) (a, F 3 , £e 3 , n) is independent of 5%. 

Now we go from H e ff >0 to H e ff (|3.3| ) considering i?ef,i as a perturbation. The twofold 
degeneracy (|3.9|) of the zero order energies forces us to use degenerate perturbation 



theory (cf. [[HJ, [28], |p9]| ). Also we have to remember that our effective Hamiltonian is 
non-hermit ian. As stated in Section T7T| , we will always assume in the general discussion 
that all non-hermitian matrices occurring are diagonalizable (For our specific numerical 
examples this assumption holds, as we see from the explicit calculations). Then the 
standard degenerate perturbation theory treatment for hermitian matrices of |30| , |28 



29] can be taken over for our case, the only difference being that projectors have to 



be replaced by quasiprojectors. 

Consider now the perturbation of the degenerate energy level characterized by 
(a, \F 3 \,£e 3 ,n) = a, 

E^ = E^(a,\F 3 l£e 3 ,n) = E w (a,-\F 3 \,£e 3 ,n). (3.10) 

The corresponding right eigenvectors of H e ff$ 

|a,±)EE|a,±|F 3 |,£e3,n) (3.11) 

span a space TZf. According to the methods of Kato and Bloch (cf. |30fl-p9|) the 
perturbed eigenvalues can be obtained by solving an eigenvalue problem in 1Z ( & . Using 
first Kato's method (cf. we define the following matrices: 

pt°> = J2 \a,m){aPm\, (3.12) 
Qg' = 1-Pg>, (3-13) 

(T)(0) P ( _ 0) 

- V 13 (k - 1 9 V 

€ fea ^ ~ ) 

(D (0) 

Sa,k = for k = 1,2, .... (3.14) 

Here the P?',^ 7^ a, denote the quasiprojectors on all the other degenerate eigen- 



subspaces of H e ff >0 , corresponding to the eigenvalues E? . We define furthermore the 
operators: 

[k%, k n+ i) := Sa,k 1 H e ff i iSa,k 2 ' ' ' Hzff^Sa^+n (3.15) 

where the kj are non-negative integers. We note the relation 

(kl, k n , k n+ i) ■ (/1, I2, lm+l) = [1 — 5o,fc„+i — ^0,ii](^l) •••> k n , k n+ \ + lx, I2, lm+l)- 

(3.16) 
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With the help of the operators ( |3.15|) we define two further matrices over TZ& as 
follows: 

oo 

k a = pl o) + £ 

n=2 

oo 

H a = E^k a + ^B' in \ (3.17) 



71=1 



where 



— '(n) 

B ' in) = E f JO,^, 0), (3.18) 

'In— 1) 



and 



^(O,^,...,^!^) (3.19) 
denotes the sum over all non-negative integers ki, k n _i satisfying 

fci + ... + k n - 1 = p. (3.20) 

In Kato's method the eigenvalues E & of H e jj = H e ff Q + H e g i are obtained as solutions 
of the following eigenvalue problem in : 

(H a - E^ka) | ) = 0, (3.21) 

where | ) G Tlf . 

We show now that both, K a and H a are matrices over which are multiples of 
the 2x2 unit matrix there. Indeed, let us first consider the matrix elements of K a : 

K a \a,m) = \a,m')(K & ) m , m , 

{Ka)m' m = (a,m'\K a \a,m). (3.22) 

From T-invariance we get, using ( |A.2j ), ( |A.3| ) and ( |A.12j ) of Appendix [A| and that n (R) 
in (|3.6|) is a real function: 

(P^h,...,^^)!) 7 = (0, ,^,0), (3.23) 

[T}k- a X) T = K a , (3.24) 

(a, m'\k a \a, m) = (a,m'\T(T^ K a T)T^\a,m) T 

= (_ 1 )2F+(m+ m ')|F 3 |(a ~m\K a \a, -m') 

(to, to' G {+1,-1}). (3.25) 

Note that the matrix T acts on the internal atomic state vectors only, but the trans- 
position has to be taken including the cm. part of the state vectors. In ( |3.25| ) F and 
| ^3 1 = —F, ...,F are half-integers. Therefore we get: 

(a, +\K a \a, +) = (a, — \K a \a, — ), 

(^+\K a \a,-) = 0, (3.26) 
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i.e. K a is proportional to the unit matrix in 7Z^ . The proof that also H a is a multiple 
of the unit matrix in 1Z { ^ is analogous. As a result we see that the twofold degeneracy 
of the eigenvalue (|3.10|) is maintained to all orders in the perturbation theory. 

The actual calculation of the perturbed eigenvalues is best done with the help of 



Bloch's method |29| which gives a smaller number of terms which have to be calculated 
in each order. Let us define the matrix H a by 



H a K a 



H a , 

Hfy. 



(3.27) 



One finds the following expansion for H a : 

oo 

where 



(3.28) 



n=l 



(3.29) 



C' w = (0,0), 

C ' (n) = E' ( n (0, A*,..., *„_!,()), (n>2). 
Here J2{ n -i) denotes the sum as in (|3.19|) , (|3.20|) , with the additional requirement 

k\ + ... + k m > m for m = 1, n — 1. (3.30) 
Clearly, the eigenvalue problem ( |3.21| ) is equivalent to 

{H & -E*)\ ) = 0, (3.31) 



with | ) G 7Z a . Furthermore H a is a multiple of the unit matrix in lZ a , as we see from 
( fT27| ). Thus we find 



E R 



|Tr H a = + ^ 



n=l 



E f = E' |Tr(0,A: 1 ,...,fc n -i,0), (n = 1, 2, . . .)• 



(3.32) 



The first three correction terms read: 



|Tr (0,0), 
Ej?) = |Tr (0,1,0), 

= iTr {(0,1, 1,0) + (0,2, 0,0)}. 



(3.33) 



Using ( 13.16| ) and arguments analogous to (|3.22|) ff. we can write E^ as follows: 



Ef = |Tr (0, 1, 1, 0) - [|Tr (0, 2, 0)1 \±Tt (0, 0) 



(3.34) 
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Now we consider a reflection R on the 1-3-plane R 2 = A 2 /2 (cf. Fig. 1): 





fx, ) 




( *i \ 




x 2 




A 2 -x 2 








K *3 ) 



(3.35) 



The reflected field £^ has its 1- and 3-components unchanged, the 2-component 
changes sign (cf. ( |1.70| )). According to Sect. [LJ] we want to calculate the following 
energy difference (cf. ( |1.25|) with £ P — > £ R ) 



En(£) — En,(£r 



(3.36) 



which is a parity-odd observable. From ( |3.4| ), ( |3.5| ) we see that for our neutral atom 
(q = 0) we have 



Heff,o(£n, Si) = Heff t o(£,6i), 



(3.37) 



H effjl (£' R ) = nH effil (£')n^. 



(3.38) 



Here 1Z is the operator representing the reflection R as defined in ( |1.72| ). From the 
R-transformation of the states ( |3.6|) as given in (|A.13|) we get the following: 



(0, ki, & n _i,0)| £R)(5l 
&a(£n,Si) 



11(0, k h k n -i, 0)\s,-6iT^, 

nk-j£,-b x )ii\ 



(3.39) 
(3.40) 



±Tr 



2 J - L 



H & (£,S 1 )-H & (£ K ,S 1 ) 



(3.41) 



Let us now have a closer look at the traces occurring in ( |3.32| ). We define the 
spherical components of the dipole operator D and the fields £ (CT) as usual: 



D 



T^(£i±i£ 2 ) 



(3.42) 



and analogously E'^ (m — ±, 0; a — 1, 2, 3). Under the reflection R (|3.35|) we have 

n*D m K = (-l) m D_ m . 

From ( |3.15| ) and (|3.38|) we get 



|Tr (0, fci, ...,k n -i,0) 

(D^ l \k 1 ,D^,...,k n ^,D^) s ■ (-1) 



(3.43) 
(3.44) 



<T,m 



mi+...+m n c ( CT i) 

^— mi 



■8. 



'(°"n) 
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where we set 



(p%\ h, D*g, h, k n ^, D^) Si (3.45) 
= (-l) n ±Tr {S^D^R^ - Rj) B(R 1 - R^- 1} ) 

x Sa,k! ■ ■ ■ Sa,k n - 1 D mn <d(R ( { n) — Ri) <d(Ri — R { { n 1) )Sa,o]- 

The point is now to realize that these expressions are invariant under rotations 
around the 3-axis performed in the internal space only. This is clear from the factor- 
ization of the states ( |3.6| ) into the internal and the cm. motion parts and the definitions 
( |3.14j ), ( p.45| ) and leads to the following results: 



5i 







for ^ rrij ^ 0. 



The time reversal transformation T gives (cf. (3.38) of [|T5| and (|3.23|) ) 



(3.46) 



(3.47) 



For n = 2 we get a further relation from T. Due to the factorization of the states ( |3.6|) 
into the internal and cm. parts we have 



(3.48) 



as is proved easily using (3.38) and (3.47) of |L5] and (|3.14j ), ( |3.4 
The reflection ( |3.35D gives 



n( CT i) u n(o-n) 



I (°0 
Rra 



Si 



(-1)"^, 

En 



(3.49) 



^ —mi i /v l) ■■■) ""n—\i - Ly — m r 



-Si 



(3.50) 



With these relations we get for ( |3.36| ), ( |3.41|) : 

oo 

AE a = E A 4 n) , 

n=l 



-Si 



X 



En I 
j=1 ™ 3 ^Vi) _ 



El 



Till 



C (.an) _ C (o-j.) C (<r n ) 



(3.51) 



where only the terms with Hrrij = are non-zero (cf. ( p. 46 )). In ( |3.51| ) the energy 



shifts AE^ are written in a way to exhibit explicitly that they are odd functions of 
Si and proportional to parity-odd combinations of the electric field components £' " . 
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Explicitly we get from ( |3.51| ) for the lowest n-values: 

|Tr(0,0) 5l -iTr(0,0)_ 5l 



El 
o. 



Si 



DP 



-Si 



C ( CT ) C 

c ~ c 



(3.52) 



AE^ = |Tr (0, 1, 0)^ - |Tr (0, 1, 0)_ 5l 

= IE[(K CT1 M^- 2, )) 5 -{d™,i,d™ 



-Si 



d { ^\ i, dX 2) ) s + (D ( y, i, d { ; 



-Si 



(°"l)|7 ( CT 2) _ £ £ 



£^l)£ 



where we used (|3.48 ). In a similar way we find 

Tr(0,2,0) 5l -Tr(0,2,0)_ 5l =0. 
The first non- vanishing result is the term with n = 3. We get from ( |3.34|) 



(3.53) 



(3.54) 



AE£> = ±Tr(0,l,l,0) 5l -±Tr(0,l,l,0)_ 5l 

- |Tr (0, 2, O^fTr (0, 0) Sl + ±Tr (0, 2, 0)_ 5l |Tr (0, 0)_^. (3.55) 



From ( |3.52| ) and ( |3.54| ) we find that the 3rd and 4th terms on the r.h.s. cancel. Thus 
we finally obtain: 



El 



n(<n) i n( CT 2) i n( CT ; 



3) 



m 3 JJ Sl 



n( CT i) i n( CT 2> i n( CT 3> 



i \m 1 +m 2 +m 3 e'Ol) p'{"l) p'(<*z) _ cVi) c'(o- 2 ) c\o$) 



(3.56) 



To show that this term is actually non-zero, let us consider the example where the 
perturbation is absent in Section o = 2: 



7/(2) 



0. 



We then get from ( |3.56j ): 

A£f = 8 1 [c«e 3 • S' w - c< 3 >e 3 • £ m ] ± (e 3 • (£' w x £' (3) )) , 
where we set for a = 1, 3: 

8 x c%> = -2i£ 3 -^(D^,l,D^,l,D^) h + (D^,l,D^,l,D^ 
+ (d?, i,D^) h - fa - -tfj)} , 
(r = 3 for a = 1; r = 1 for a = 3). 



(3.57) 



(3.58) 



Si 



(3.59) 
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To give some numbers we consider a box (Fig. 1) of dimension 



A x = A 2 = A 3 = 8nm (3.60) 

and segment lengths (cf. Q2.2|) ): 



R[ 1} = 0.4 • A u i?f - 0.6 • Ai. (3.61) 

In Fig. 2a we plot the real parts of the corresponding unperturbed eigenenergies 
( p.lO| ) for our "hydrogen" atom with I = as function of the electric field £. We note 
that the transverse parts of the wave functions, characterized by the integers ri2,n 3 
(cf. ( |3.6| )) do not mix under perturbations by H e ff t \ ( |3.5| ) for our case here (q = 0). 
Thus we consider only the levels with the lowest transverse energies corresponding to 
ri2 = 11% = I. For the Lamb shift L and the 2P3/2 — 2P1/2 energy difference A as well 



as the widths T$, Tp we take the values from real hydrogen (cf. Table III of JT5|). The 
zero of the energy scale is put at the energy of the 2Pi/ 2 level in vacuo. 
Now we introduce the perturbation in the electric field 





(3.62) 



This leads to an energy shift — > E& as calculated in ( |3.32| ) and to a P-violating 
energy difference AEa- This is given to lowest nontrivial order in ( p.58|) , from which we 



see that AE& } is proportional to r)S' 3 and thus changes sign if we change the chirality 
of the field arrangement (77 — > —77): 




A££>=<f lC gM- • ( 3 - 63 ) 



We plot the real and imaginary parts of the coefficient Cq for the state a = 
(25i/2, 1 F 3 1 = l/2,ni = 1,71,2 = 1,713 = 1) in Fig. 3. Of course, we would like to 
make |Ai?^ 3) | as big as possible by increasing £'. However, there are limits to this 
in our calculation. The perturbation treatment is only reliable as long as the energy 
shifts induced by the perturbation are small compared to the spacing of the unper- 
turbed levels. Thus, we can at best trust our calculation for values of \£'\ < £' max 
which we estimate from the requirement that the first plus second order energy shift 
is equal to the minimum of the unperturbed energy differences: 

W(±CJI + W(±C»)I = mj* K ~ Eg>\. (3.64) 

Here we choose the plus or minus sign according to which gives the lower value for £' max . 
We plot £' max in Fig. 4 for the level a = (2S 1/2 , \F 3 \ = 1/2, m = 1, n 2 = 1, n 3 = 1). The 
P-violating energy and lifetime shifts for £' = £' max are then shown in Fig. 5. From this 
calculation we can thus establish theoretically for our "hydrogen" atom with / = a 
P-violating energy shift of the order of 

^AEf/h ~ • 5 • 10 3 Hz = -2.46 • 10~ 9 Hz. (3.65) 
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Since Si is of the order 10- n ...10- 12 (cf. (3.20), (3.21) of fT§), this is an extremely 
small energy shift. Fortunately there is — at least in theory — a way to enhance the 
effect dramatically. In fact, we will show now how to obtain energy shifts not of order 
Si but of order y/S~i. 

Consider Figure 2a. There is a crossing between the levels a = (2Si/ 2 , \F 3 \ = 
1/2; n = (1,1,1)) and (3 = (2P 1/2 , \F 3 \ = 1/2; n = (2,1,1)) in the real parts at an 
electric field strength £ ~ 500 V/cm. Correspondingly, the maximally allowed pertur- 
bation field strength £' max is at a minimum at £ ~ 500 V/cm as can be seen in Fig. 4. 
This can be easily understood from ( |3.64| ) and the fact that l^E^ /h\ ~ 16 MHz and 



\^E^/h\ ~ 38 MHz are by two orders of magnitude smaller than \^tE ( °^/h\ ~ 1.6 GHz 



(cf. Figs. 2a, 2b). Then, the r.h.s. of Q3.64 ) becomes very small at the zero of the real 



part ^(-Ei 0) — E^) and so does £' max . However, this only means that our calculation 
as done so far is limited. In fact, the P- violating energy and lifetime shifts at such a 
crossing of levels may nevertheless have a maximum. We will show this in the following 
using again degenerate perturbation theory treating all the crossing levels together. 

We consider a configuration with a crossing of levels in both the real and the 
imaginary parts at the same field strength: 

M(E$\£e 3 )-E^(£e 3 )) = 0, (3.66) 
Z(E$\£e 3 )-E^(£e 3 ))=0. (3.67) 

We find (cf. Fig. 2b) that ( |3T67D is fulfilled at 

£ res ~ 1024 V/cm (3.68) 

for a = 2Si/2, (3 = 2Pi/ 2 . If one now adjusts the box lengths to 

A x = 5.84 nm, A 2 = A 3 = 8 nm, (3.69) 



then ( ETE5D is satisfied for a = (2Sy 2 , \F 3 \ = 1/2; n x = 1),(3 = (2P 1/2 , |F 3 | = 1/2; m 



2) and so is (|3.67|) as there are no contributions to the imaginary part of the energy 
from the cm. motion. Note that ( [3.661 ) could also be fulfilled by adjusting A 2 , A 3 , but 



the differences between energy levels with different cm. quantum numbers n 2 , n 3 will 
not appear for neutral atoms as energy denominators in the perturbation series, since 
the perturbation fields £' {a) are chosen not to vary along the 2- and 3-directions. Then, 
due to orthogonality, only matrix elements 

(a,F 3 ; (ni,n 2 ,n 3 )\H eff> i\P,F^ (n[,n 2 ,n 3 )) 

diagonal in n 2 ,n 3 are nonzero. 

In order to calculate a reliable value for the P-violating energy shifts, we use de- 
generate perturbation theory in the 4-dimensional subspace spanned by the states 

|«,±) := \2S± ±|F 3 |,£ res e 3 ; (1,1,1)), |A±) := \2P± ±|F 3 |, £ res e 3 ; (2, 1, 1)), (3.70) 
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where \Fs\ = 1/2. In analogy to ( |3.12| )-( |3.14| ) we define the matrices: 

\a,m)(a,m\ + \/3,m)(/3,m\ 



p(0)_ 

© (0) - 
© (0) - 



s, 



(a/3),0 



S, 



(a/3),k 



E 



(a/3)' 



E 

_p(0)_ 

(D (0) - 

for fc = l,2 l ... 



(A; = 1,2,...) 



(3.71) 
(3.72) 

(3.73) 



(3.74) 



(Remember that in ( |3.73| ) = E^ ^ E^ ). Using these matrices we define, similarly 
to ( |3.15|) , the operators 



(ki, k n +i) :— S^p),k 1 H e ff,iS t 



(&p),k 2 ' 



n+l 



(3.75) 



In the following we will show that all energy levels remain at least twofold degenerate 
to all orders of perturbative theory also in this case. We define 



(a/3) 



(a/3) 



/(«) 



n=2 



H, 



(a/3) 



(3.76) 
(3.77) 



n=l 



with v4 /(n) , B' {n) defined as in ( 3.18Q , but in terms of the operators (|3.75|) . 

Instead of discussing now Kato's eigenvalue problem ( |3.21| ) we will now introduce 
and use a modified version of Bloch's method |29|. We define a matrix H^p) by the 
conditions (see Appendix |C|): 



TT _ ip(Q) 
]P(0) TT _ 



H, 



(a/3)> 



H, 



(a/3> 



(3.78) 



Bloch's operator Ht&m defined as in ( p.27| ) with a —>■ (a/3) is in general non-hermitian 

even if H( & p) an d K(ap\ are hermitian. In contrast, H^p) is hermitian if H<up\ and K^p) 
are hermitian, which would be the case for stable states, i.e. T$ = IV = 0. We will 
need this property of H^p) later when we want to check, that our effect proportional to 
\f5\ vanishes for the case of stable states. We assume that the perturbation expansion 
of H(ap~) is well defined and converges (cf. Appendix for the conditions for this to be 
the case). The eigenvalue problem Q3.21| ) is then equivalent to the following standard 
diagonalization problem: 



H, 



(a/3) 



E)\ ) = 0, 



(3.79) 
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where I ) G 

I I (a/3) 

Now we will show, that the eigenvalue problem factorizes and leaves a twofold 
degeneracy of the eigenenergies. We define the 2x2 matrices: 



([i,m'\H {dl p)\l,rn) 
where 7', 7 G {a, /?}; m', m G {+1, — 1}. Now the secular matrix reads: 



S(E) := ^ Ef 

By simple manipulations we obtain that 

detS{E) =det S\E) 



E 



E 



where 



S'(E) 



(3.80) 

(3.81) 

(3.82) 
(3.83) 



H a pHp a H a p(Hpp — E)Hp & 

Using T-invariance we show in Appendix that Mm and Ma^Mpa are all mul 
tiples of the 2x2 unit matrix. 

Then it follows from (|3^ ), that 



det S(E) = ((a - E)(b - E) - c 



(3.84) 



where 



a 
b 

c 



-Tr H-- 

2 — aai 

|Tr Hpp, 

^Tr {HafiHpa). 



(3.85) 



Thus det S(E) is a perfect square, all eigenvalues E are twofold degenerate and 
det S(E) = if and only if 



(a -E){b-E) -c = 0. 
This is a quadratic equation for i£ which has the solutions 



(3.86) 



E 



■■ \ a + b± ^{a-b) 2 + Ac 
= \ a + b±yj[-i(a - b) + 2y/c][i(a - b) + 2^] 

The complex numbers a, 6, c can be expanded in 5\\ 

a = do + Sidi + ... 
b = b + Sibi + ... 
c = c + 5iCi + ... 



(3.87) 



(3.88) 
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and thus we get for the eigenenergies ( |3.87| ): 



E, 



\ [oo + b + 5i{a x + 61) + 0(5f) 

±{[-i(a - b ) + 2 v /c^][i(a - b ) + 2y/c^\ 
+o 1 [2(a - 6 Q )(a 1 - b x ) + 4cJ + 0(^)} 1/2 
Now it is clear that in the case 

(a -6 ) 2 + 4c ^0 
we will get at best P-violating energy shifts linear in S\. If, however, 

i(a -b ) + 2y/c^= 



or 



- i(a - b ) + 2^ = 
and together with ( |3.91| ) or ( |3.92| ): 

2(a -b )(a 1 -b 1 )+Ac 1 ^0, 



(3.89) 



(3.90) 



(3.91) 



(3.92) 



(3.93) 



we will obtain 



E 



aP,± 



a + b ± y^{2(a - &o)( a i - & i) + 4c i} 
+0(5 1 



1/2 



(3.94) 



This is a P-violating splitting which can be distinguished from a P-conserving one as 
it changes its phase by a factor of i under a change of sign of 5i or a reversal of the 
handedness of the external field arrangement which changes the signs of a±,bi,ci. In 
other words: we have to search for a degeneracy of two complex energy levels (|3.87|) in 



an external electric field with non-zero handedness in the absence of P-violation, i.e. 
for 5\ = 0. Then the small P-violation effects proportional to 5i in the matrix elements 
can produce an energy shift of order \f5[. 

Note, that this is clearly an effect which results from the non-hermiticity of Hi a p\, 
i.e., which is only possible for unstable states. For the case of stable states, H^p) is 
hermitian (cf. ( |3.78j )) and therefore we get for the submatrices ( |3.80| ) 



H 



PP 



H 



H 



a/3 



H 



13a- 



(3.95) 



From (|3.85| ) we find then that a, 6, c and likewise ao, bo, Co, and ai, b±, c\ are real and 



c > 0, c > 0; 
ci = if Co = 0. 



(3.96) 



Now the condition (|3.91|) is equivalent to (|3.92|) and leads to ao = b and Co = which 
implies through ( |3.96D also c\ = 0. Therefore, ( |3.93|) cannot be fulfilled simultaneously 
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with ( |3.91| ) or, equivalently, Q3.92| ) for the case of stable states, and there will be no 
energy shifts proportional to \f5[ for them. 



On the contrary, if a, b and c are complex numbers, which is the case for unstable 
states, one can in general fulfill ( |3.91|) or ( |3.92|) together with (|3.93|) . We will show by 
the results of a numerical investigation, that this is indeed the case. We choose the 
following general configuration of three perturbation fields: 



(<0 



£'; 



1 o 



1,2,3). 



(3.97) 



In Appendix |D] we give a proof that for sufficiently small £' there exist indeed "reso- 
nance" values for the vectors T7 (CT) where energy shifts proportional to y/S\ occur. As 
an example we choose £' = 3 V/cm. Then the calculation to first order in £' gives the 
following resonance values: 



'Ires, 




(2) 
'Hres, 



( 








„(3) 
'Ires, 



\ 574.93 




(3.98) 



where (|3.92 ) and ( |3.93 ) are fulfilled. Higher orders in £' can only shift these values 
slightly. 



In Figs. 6 a and b we plot the decadic logarithms of the moduli of the real and 
imaginary parts of the P- violating energy shift (cf. fl3.94j )) 



W{£ V) }) - W{£r (ct) }) 



(3.99) 



for our "hydrogen" in a field configuration given by ( |3.61| ), ( |3.68j ) and ( |3.98j ) vs. the 
deviation -q^ — ri ( ^ res _ of 1]^ from the resonance value r]^ res _ = 574.93 for which 
the P-even splitting of levels is removed, revealing an enhancement of the P-violating 
splitting of the order of l/V^i- The width of the peak is proportional to \f5[ an d 
found to be A^ 2 ' ~ 0.4 • 10~ 4 , corresponding to £' Ar/^ = 1.2 • 10~ 4 V/cm. Note, 
however, that in Fig. 6a the differences of the real parts of the energies of the levels 
are at maximum of the order of 10~ 2 Hz whereas the lifetime of the levels is calculated 
to be 0.29 ■ 10~ 8 s. This corresponds to a line width AE/h = 5.44 ■ 10 7 Hz. Thus it 
is certainly not an easy task to detect frequency shifts of order 10~ 2 Hz in such broad 
lines. But our calculation is only meant to provide an illustration how energy shifts 
proportional to y/5[ can be obtained. Certainly, for a realistic experimental situation 
one would try to optimize such an effect. 



3.2 Results for an ion with 1 = 

In this section we will discuss the modifications one has to include for the physically 
relevant case of ions without nuclear spin, as 2He + . As we have already discussed in 
Section 2/2 one has to introduce Airy functions instead of sines and cosines for the cm. 
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wave functions. We write the (right) eigenstates of H e ff$ as 
|a,^3,£e 3 ,n) := \a, F 3 , £e 3 ) • n (R), 

n (R) := AA n3 ^+(^3,^3,0 -^-(^3,^3,-^)] 



2 



2 ( n^R 



xnvx sin PP)' (3 - 100) 



where the Airy functions ^ are defined as in ( |2.37] ) , ( |2.38| ), Q2.40| ) f. The eigenener- 



gies E n3 corresponding to the Airy modes, which are numbered by n 3 , are obtained by 
solving a determinant equation analogous to ( |2.43| ), the coefficients £ n3 and the normal- 
ization A/"„3 follow from (|2.40|) and Q2.44|) , respectively Using the eigenstates (|3.100|) 
one computes the perturbation series in an analogous manner as described in ( |3.12| )- 
( p.20| ) for Kato's and ( |3.27f )- (|3.30| ) for Bloch's method. Note that the perturbation 
H e ff i ( |3.5| ) to the effective Hamiltonian i?eff,o now contains the additional term 

3 

- 9 R^£' W 9« -R 1 )Q{R 1 -R^), (3.101) 

(7=1 

which was absent in the q = case. Therefore, the terms of order n = 1,2, ... of the 
series ( |3.17| ), (|3.18|) , (|3.29| ) have the following structure (symbolically) 



A' (n) oc ((qR + UT) 

oc ((D) n + n(D)"- 1 (gR) + ...) (3.102) 

and similarly for B' {n \ C /(n) . 

Since matrix elements of (|3.101|) w.r.t. the states (|3.100|) are diagonal in the internal 
quantum numbers (a, _F 3 ), ^i-linear P- violating contributions to the series can only arise 
from matrix elements of D. Therefore, up to the third order, n < 3, there will be no 
P-violating contributions to the energy shifts AE& (cf. ( |3.52| ) ff.) involving matrix 
elements of (|3.101|) . This follows from our result that such contributions appear only 
for n > 3 in the q = case. Hence, the expressions ( ^.52[ ), ( |3.53p , ( |3.58| ) f. for the 



P-violating energy shifts AE a up to n = 3 remain valid for ions with 1 = 0. 

However, when one calculates the maximally allowed perturbation field strength 
£' max using ( |3.64j ), the contributions from the matrix elements of ( |3.101| ) to the P- 



conserving energy shifts E & ,E a have to be taken into consideration. 

We will now present results for |He + -ions in a box of lenghts: 

A 1 = A 2 = A 3 = lnm. (3.103) 

Note that compared to the case of our "hydrogen" (cf. ( p.60| )) we rescaled the Ai by a 
factor 

Z- 2 {Z + N)- 1/2 = 1/8 (3.104) 

for ^He" 1 " (Z = N = 2) in order to obtain an approximate scaling factor of Z 4 for the 
cm. energies. For £' = the contribution of the c.m.-mode in i-direction reads 

E ni = (3-105) 
2 m Af 
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As the mass of the nucleus of a light atom may be approximated by (Z + N) times 
the mass of the proton, we obtain the desired scaling factor Z A for the cm. energies 
Q3.1051 ), if we multiply the widths A4 of the box used for hydrogen with (|3.104|) : 



^n 2 Z 4 [Z + N] 4 



As it was mentioned in section 3.2 of |15|, one obtains "reduced Z = 1" problems 



from the general Z case if one multiplies all electric field strengths of the Z = 1 case 
by Z 5 . In the following we will choose such rescaled fields in addition to the rescaled 
box widths and observe, that all energies will approximately be larger by a factor of 
Z A = 16 compared to the results for "hydrogen" . 

In Fig. 7 we plot the lowest levels of the cm. motion in 3-direction vs. the electric 
field strength £, which are the solutions to ( 2.43|) , the eigenenergies of the Airy states. 



As the zero of the energy scale we choose the "minimum" of the energy in 3-direction 

F - / ° for ^<0, 

E ^-\ - qAs S for £>0. {3A07) 

For £ = the quadratic dependence on n 3 corresponding to the sine-/cosine modes is 
recovered. 

The Lamb shift L and the energy difference A as well as the widths Fs, Fp for ^He 4 " 



are given in Table II of ||15|| . In the following we again restrict ourselves to levels with 
the lowest cm. energies corresponding to n« = 1, i — 1, 3. 

We plot £' max in Fig. 8 for the level a = (2S 1/2 , \F 3 \ = §,n = (1,1,1)). In Figs. 
9 a, b we show the P- violating (a) energy shifts and (b) shifts of the decay width for 
£' = £' max - We have chosen the same configuration of the perturbation fields as in 
( gtflp , ( gjjD with 77 = 1. 



From our calculation we obtain a theoretical energy shift of the order of 

ZtAEg/h ~ 5 1 ■ 1.2 ■ 10 5 Hz = 1.36 ■ 10~ 6 Hz, (3.108) 



which is larger by a factor of 550 > Z A 5i(2, 2)/#i(l, 0) compared to the result (|3.65|) 
for the "hydrogen atom" with 1 = 0. 

As one can expect there is the same possibility as in the q = case to obtain 
P-violating energy shifts of the order of y/5[. 

The calculation is the same as described in Section |3T, (|3.66|) - (|3.971) , and Appendix 



0, piD-flrrSD and ( |TXT3D f., except that it now involves the states (PIDOP and the 
additional piece ( p,101| ) to i?ef,i- The secular matrix may be block-diagonalized in 
the same way as in ( |3.8U| )— ( 3.84|) because ( |3.101| ) is diagonal in the internal quantum 



numbers (7, m). 

We have chosen the crossing of the levels a. = (2Si/ 2 , n = (1,1,1)),/? = (2Pi/ 2 , n = 
(2, 1, 1)), (jOED , (EUTTp , which occurs at: 

£ res = Z 5 ■ 30.11 kV/cm (3.109) 
for A x = 0.76 nm, A 2 = A 3 = 1 nm . (3.110) 
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For the case of a 2He + ion in a box described by ( |3.61| ), ( |3.109| ), ( |3.110| ), we choose 
again perturbing fields as in (|3.97|) . We show in Appendix |D| that for 



£' = 96 V/cm. 
the resonance values for r) (,T) to first order in £' are 



(3.111) 



'Ires 



1 



-47.37 



? "Hre 



(2) 






356.97 



„C3) 
5 'Ires, 





1 

-49.73 



(3.112) 



Here again the analogues of (|3.92|) and (|3.93|) are satisfied. 



In Fig. 10a and b we plot the decadic logarithms of the moduli of the real and twice 
the imaginary parts of the P-violating energy shift 

E^ + ({£'^}) - E^ + ({£'^}) (3.113) 

for pie + in a field configuration given by ( |3.109|) and ( |3.112|) f. vs. the deviation rj 3 — 
V3,res,- °f Vtt from the resonance value ry£ res _ = 356.97 for which the P-even splitting 
of levels is removed, revealing an enhancement of the P-violating splitting of the order 
of l/^/6l. The width of the peak is proportional to y/5[ and found to be A?7 3 2) = 1CT 3 
(S'Ar]^ = 9.6- 10~ 2 V/cm). In Fig. 10 a the differences of the real parts of the energies 
of the levels are now at maximum of the order of 1 Hz. The lifetime of the levels is 
calculated to be 0.19 • 10~ 9 s, which corresponds to a line width AE/h = 8.60 ■ 10 8 Hz. 
As before it will not be easy to detect frequency shifts of order 1 Hz in such broad lines. 
Also this calculation is only meant to show that energy shifts proportional to y/Si can 
in principle be obtained also for ions and it illustrates how q 7^ affects the numerical 
solutions. 



3.3 Results for a neutral atom with 7 = 1/2 

Let us now turn to the case of a neutral atom with / = 1/2, i.e. JH. As discussed in 
Appendix [A], the eigenstates of H e ff iQ (£e 3 ) ( |3.4j ) are twofold degenerate in F 3 = ±|F 3 | 
for fixed l,j,F, if F 3 ^ 0: 

E m (a,F 3 ,£e 3 ,n) = E m (a, -F 3 ,£e 3 ,n), (3.114) 

where a = (2, l,j, F), but there is no degeneracy between the levels of different F, 
but equal \F 3 \. Hence, for the calculation of perturbative corrections to levels with 
I-F3I 7^ 0, the same method involving degenerate perturbation theory applies as it is 
described for the 1 = case in Section |3.1| . As far as the | ^3 1 =0 levels are concerned, 
simple non-degenerate perturbation theory will be sufficient. 

Let us concentrate first on the question, if the F 3 degeneracy of the |F 3 | 7^ states is 
in general broken by the perturbation fields £ l{a) . Our results of Chapter [2.1.1| indicate 
that this will be the case. 
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Consider the degenerate energy levels characterized by a = (a, \F^\,Se^, n), | ^3 1 7^ 
0. The corresponding right eigenstates of H e ff t0 (( |3.4| ) with q = 0), 



\a, ±) := \a, ±\F 3 \,£e 3 , n) 



(3.115) 



span a two-dimensional space IZ^ . Now we turn to #eff,i (Q3.5D with q = 0) and use 
the same notation for the perturbation theory treatment as in ( p.l2|) -( [3.21| ). 

From T-invariance we obtain again the relation ( |3.25| ). However, now F and | ^3 1 
are integers and we get 



(a, +\K a \a, +) 



at, 



\K a \a, -), 



(3.116) 



and self-identities for the off-diagonal elements, which therefore need not to be zero. 
The analogous relations apply to the matrix elements of H a . 

For both, the general analysis and for actual calculations we use our modified Bloch 
method defining H a in analogy to (|3.78|) (cf. Appendix |C|) . Then, in order to solve the 
eigenvalue problem analogous to ( |3.79| ), one has to solve the secular equation 



det (H 6 



0. 



where 



a, m' \ H n \a, m 



)) - {H. 



From T-invariance we find, similarly to (|3.116|) : 
The solutions of ( |3.117| ) are: 



E a ,± — H a++ ± 



H a ,^ — ELa, — h 



1/2 



(3.117) 
(3.118) 

(3.119) 
(3.120) 



Thus, in general the twofold degeneracy of the levels of -ff e #,o f° r I ^3 1 7^ is broken by 
the perturbation. 

Let us now consider the perturbation expansion of H a according to ( |C.8| ) of Ap- 
pendix |C[ 



E a = E^ + EC 



(n) 



(3.121) 



n=l 



For the diagonal elements H a mm the same arguments apply which are given in Section 

3~T| (( p.42|) -( p.56|) ) for the trace of H a . In this way we see that there will be no ^i^-linear 
contributions to H a ^ mm below the third order. 

As far as the off-diagonal elements H ajm i m , m! 7^ m are concerned, they are non-zero 
only above the first order. Indeed, for the first order term in (|3.121|) we have: 



(a, m'\ C {1) I at, m) = (a, m'\ i?eff,i | or, m). 



(3.122) 



For m' 7^ m this matrix element vanishes due to the same argument which led to (|3.46|) 
and the fact that the difference AF3 in F3 between the degenerate unperturbed states 
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( |3.115| ) satisfies IAF3I = \m — m' 1 1 1 > 2 and that the dipole operator ( |3.42| ) cannot 
change F3 by more than one unit. 

We consider now in detail the following configuration of perturbing fields: 



E' w : = rj M • £', (a = 1,2,3), 



(3.123) 



where the dimensionless vectors 77 (ct) are left arbitrary for the moment. For the expan- 
sion of Ha t m'm we have then from (|3.121|) : 



with 



E^S m , m +J2 («,m'|C ( ' l) |«,m), 



n=l 



(3.124) 



C {n) \a,m) 



From the discussion above we find 



KO) , > . > 
a,m'm 1 / - ' - a,m'm 
\ i=l 



m',m G {+1, — l},n > 1. 



„(»>*) 

„(M) 



for n < 2; i = 1, 2, if m' = m; 
for z = 0, 1, 2, if m ^ m. 



(3.125) 



(3.126) 



Using these relations we can easily convince ourselves that a consistent evaluation of 
E & ,± Q3.120Q up to order E' 2 requires only the knowledge of the coefficients c^/ m for 
n < 2. We get then: 



n=l 

/ r 2 

(2,0) , V^r (2,0 



± 



i=l 



,(2,0) 



£ ,2 + 0(£ ,3 ).(3.127) 



In the general case where all coefficients 4™^/ m appearing in (|3.127|) differ from zero, 
one therefore gets for the level shifts: 



E, 



Ek 0) + E^ + R 



(2,0) 



± \ (aE^ + £ SiAEi A + 0(51^ + 0(E' 3 ). (3.128) 

In E& ,1 \ AE^' z) n denotes the order in £' , i = the order zero in 5\ y 2 and i = 1, 2 the 
order one in <5 l5 5 2 respectively. We see that here linear terms in 5 12 occur already at 
order E' 2 . 

To give some numbers we choose an electric field configuration (|3.123|) with 



(ct) 

m 



0; for cr=l,3 and rj 



(2) 



0. 



(3.129) 
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Then the first order level shift vanishes, E, 



(1,0) 



and we have 



Ea,± = ^ 0) +fE{Ci« ) < ) +^ ) ^ ) )} 

\ cr,r 

E{c 2 pK (^V T) -^ 2 T> )(<V 

CT.T I 



2 



1=1 



5 



/2 



+0(£' 3 ), 



(3.130) 



where the 5^,5^ are complex numbers which depend on S. The sums over a,r,p,K 
run over the segment numbers 1,3. In order to compute the P-violating energy shift for 
the levels a = (a,\F 3 \ = l,n = (1,1,1)) for a = nljF e {2S±1, 2P§1, 2P§1, 2P§2}, 
we take the values for L, A, Ts, Tp and the hyperfine splittings A±, A 2 , A% from Table 
III of [15] and chose the following configuration:^ 



A 1 = A 2 = A 3 = 8 nm, 
RV/At = 0.45, Rf ) /A 1 = 0.55. 



(3.131) 
(3.132) 



As one can see from ( 3.130| ), the * 12 Tinear contributions vanish if rf[' = or rj 2 ' 



(°-) 



0; cr = 1, 3; or if rji } = rj 2 CT) = 0, a — 1, 3 and 772 = ±r/j 4 CT) . We choose 



77 



ID 




77 



(2) 



0, T) 



(3) 




(3.133) 



in order to maximize the P-violating effect. In Table I we give numbers for the P- 
conserving as well as the P-violating contributions to E& y ± for different levels a at 
different electric fields S. From the analogue of ([3.64 ), we found the maximally allowed 
perturbation field strengths and used these to compute the results in Table I. 
We establish thus theoretically a P-violating energy shift of the order of 

^(E(S,{£' ia) })-E(S,{£' R (a) }))/h = ( 1.2 <f x - 20.2& + 0([S'/S} 3 )) kHz 

= -2.56 • 10" 8 Hz, 

23 (E(S,{£ ,ia) }) - E(S,{e' R {a) })) fh = (-I.U1 + 16.3 5 2 + 0([S' /S} 3 )) 10V 1 

= 2.08 ■ 10~ 7 s" 1 (3.134) 



for a = 25*1 1 where the numerical values for * 12 are taken from Table III of |Tj|. We 
observe that the nuclear spin dependent term oc 5 2 of the energy shift is enhanced 
by a factor of ~ 17 w.r.t. the spin independent one. As was already noted in 



2 In this paper we use the notation Ai, i = 1,2,3 for the cartesian lengths of the box, which are 
not to be confused with the hfs-splittings which we denote here by A$, i — 1,2,3. 
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an enhancement of the terms proportional to 62 over the ones proportional to 81 is 
expected to be a general feature of P- violating effects for {H. (cf. (6.39), (6.40) of [|1^ 



Numerically an energy shift (|3.134|) is, of course, very small and presumably not easily 

measured. 

Let us now consider the possibility of P-violating energy shifts of the order of \J$i,2- 

From (|3.120[) we can easily see how this can be achieved. Let us expand Hu,m'm m 6\^- 



Br, 



H, 



(()) 



(1) 

a,m'm 



+ S 2 K 



(2) 

a,m'm 



(3.135) 



For unstable states as we are considering here H a is a non-hermitian matrix and we 
can have 



rr(0) 
n a, -f 

rf(i) 

H { - 0) 



= 0, 
* 0, 
± 0. 



1,2, 



(3.136) 
(3.137) 
(3.138) 



In this situation the level shifts ( |3.12U| ) read: 



E, 



a,± 



rr(0) 
n a, ++ 



± 



C(5i,5 2 ). 



H. 



(0) 



1/2 



(3.139) 



Without P-violation, i.e. for 61 = 62 = 0, the levels would be degenerate. 

Note that we need not to have a crossing of the energies and lifetimes of the un- 
perturbed states in order to achieve a degeneracy of levels for nonzero handedness of 
the external fields. As the degeneracy of the ±1^1 levels gets broken by electric field 
perturbations anyway, it suffices to suppress this breaking by a convenient choice of 
77 (ct) , a = 1,2,3 at any field S. Hence, we can choose a convenient zero order field 
strength of e.g. 

S = 200 V/cm. (3.140) 
The condition (|3~T36D reads now (cf. fl312l) - (|3TT26D ): 



H, 



(0) 



n=2 

(2,0) C-/2 1 (3,0) C-/3 , 
Lq -| — O T L-ci -| — O ~r • • • 

0. 



To order £^ /2 we can easily solve this equation. We found a solution to 

r (2,0) _ Q 

to be 

/ 3.91 \ / 3.91 



77 



(1) 



0.0764 I , T? (2) = 0, r/ 3) 

V 



-0.328 




(3.141) 
(3.142) 

(3.143) 
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where we chose the box geometry as in ( |3~T3"T| ), ( gTgg) . Note that to order £' 2 the 



magic condition ( |3.142|) does not depend on £' and we get P-violating energy shifts pro- 



portional to ySi,2 for all £'. However, the requirement of applicability of perturbation 
theory (cf. (|3.64 )) gives here as maximally allowed value 



Ca, = 8-3 V/cm. (3.144) 



The reason for this relatively small number is that the r.h.s. of ( 3.64|) is now of the order 



of the hfs splitting, which, at £ = 200 V/cm, is by two orders of magnitude smaller 
than the differences between levels for nuclear spin 1 = 0. Thus, choosing for instance 
£' = 5 V/cm < £' max we are safely in the region where we can apply perturbation 
theory. 

Using now the implicit function theorem in a way analogous to ( D.7 ) ff. of Appendix 



[B| we establish that a solution to ( |3.141| ) exists to all orders in £' for sufficiently small 



£'. In other words: For such small enough field strength £' higher order terms in the 
perturbation expansion ( |3.141| ) will only displace the "resonance" position slightly but 
cannot change the qualitative features. Thus we establish theoretically the existence 
of P-violating energy shifts proportional to \J$i,2- 

The signal to look for in an experiment is the difference of the eigenenergies E &t ± 
( 3-139| ) for the perturbation fields £' (a) and the reflected ones £ R <CT) 



£' R M = t,<? ■ £', (3.145) 



where rj^ are the same vectors as in (3-143 ) but with the opposite sign of the 2- 
components. 

In Figs. 11 a, b we plot the real and imaginary parts of these energy differences 

%({£' W })-%({5' R W }) (3.146) 

for £' = 5 V/cm versus the deviation of £ 2 (3> from the resonance value £ 2 ( rL = —1-64 
V/cm. Numbers for other values of £' are easily obtained from Figs. 11a, b, since (in 
our approximation) these energy differences are proportional to £' 2 . The width of the 
peak is proportional to \f~5~i and found to be 2.6 • 10 -12 V/cm. The differences of the 
real parts of the energies shown in Figs. 11a, b are at maximum of the order of 10 -5 Hz 
whereas the lifetime of the levels is found to be 4.13 • 10~ 8 s. This corresponds to a line 
width AE/h = 3.86 • 10 6 Hz. Again we note that detecting frequency shifts of order 
10~ 5 Hz in such broad lines is certainly not easy. However, our example shows, how 
energy shifts proportional to \J$i,2 can be obtained for atoms with I = 1/2. In Figs. 
11a, b we observe the rise of the signal by a factor ~ l/v^ ~ 10 6 at the resonance 
value. We are sure that for a realistic experimental situation one can optimize such 
effects. But it remains to be seen what values of AE one can then obtain. 
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4 Conclusions 



In this article we have investigated parity violation in hydrogen- like atoms due to Z- 
boson exchange between the lepton and the quarks in the nucleus. In particular we 
have considered the consequences of this P-violating neutral current interaction for the 
eigenenergies of the atoms which are moving in inhomogeneous external electric and 
magnetic fields in finite regions of space. 

In this first part we have discussed the general formalism to treat a non-relativistic 
atom in an external field arrangement. We derived the consequences of T-invariance for 
the eigenenergies of the atoms. In general these eigenenergies arise from an interplay of 
the internal and the cm. motion of the atom. We have shown that for an atom confined 
to a finite region of space where we have a homogeneous electric field the eigenenergies 
get no P-violating contribution. Then we studied configurations of an external electric 
field within a rectangular box, where the field is homogeneous throughout different seg- 
ments of the box and shows sudden variations at the interfaces of these segments. In 
part II we will study configurations where the field varies adiabatically along a certain 
direction in space. 



We studied hydrogen like atoms in states of principal quantum number n = 2 in 
detail. First we considered boxes with two segments only, where the two electric field 
vectors together with the vector which is normal to the interface form a chirally non- 
symmetric configuration which could in principle be sufficient to exhibit P-violating 
effects. However, we showed, that for this geometry there are no P-violating energy 
shifts due to the T-symmetry at least in our scheme of approximations. Thus we found 
that a box with a minimum of three segments is needed if the electric field varies along 
one direction of the rectangular box only. 

We studied then explicitly hydrogen-like atoms with vanishing nuclear spin 1 = 
like %Re + and with nuclear spin 1 = 1/2 like JH in a box with three segments us- 
ing the methods of Rayleigh Schrodinger perturbation theory. We indeed found P- 
violating contributions to the eigenenergies of the systems for field configurations car- 
rying nonzero chirality. This opens the possibility of determining the parameters of 
the neutral current P-violation in atoms by frequency measurements. 

We found energy shifts linear in the P-violation parameters 5i (i = 1,2), where 
5{ ~ 1CT 12 are the mixing amplitudes between S and P levels due to the nuclear spin 
independent and dependent P-violating contributions to the Hamiltonian divided by 
the Lamb shift determining the order of magnitude of P-odd effects. Such energy shifts 
were found to be of the order of 10 _5 ...10~ 9 Hz for boxes of ~ 5-10 nm size. To measure 
these extremely small effects would require an accuracy which only in the context of 
the hydrogen 21 cm-line-measurements, cf. 0, has been reached so far. 



There is some similarity of the energy shifts linear in the P-violation parameters 
for the systems we considered with the P-violating energy difference between the two 
enantiomers of chiral molecules. Such energy differences are studied theoretically e.g. 
in |J4, 35 1 . But we emphasize that our physical systems are very different from chiral 
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molecules and also the physics of our P-odd energy shifts is quite different from that 
in chiral molecules. In chiral molecules the P-odd energy difference occurs with no ex- 
ternal field due to a P- violating distortion of the electronic wave functions in the chiral 
arrangement of the nuclei. The cm. motion of the molecule is irrelevant for this. In our 
case we consider simple atomic systems where questions of "spontaneous" occurence of 
handed molecules, maybe due to a spontaneous symmetry breaking do not occur. 
We get P-odd energy shifts only in the presence of external fields. The interplay of the 
cm. and internal motions of the atom is crucial. Since we have only one nucleus and 
one electron around it our energy shifts certainly do not occur due to the movement 
of the electron in a chiral lattice of nuclei. As will be seen even more clearly in part 
II of this work our energy shifts arise in essence through P-violating rotations of the 
atom as a whole in the external field (cf. also [15|]). Through the boundary conditions 
for stationary states such rotations lead to P-odd energy shifts. However, we can make 
a somewhat loose connection between P-odd energy shifts in chiral molecules and our 
results in the following way. We start from our situation of an atom in a box with 
a chiral electric field and — in a "Gedankenexperiment" — imagine to make the box 
smaller and smaller until it has atomic dimensions, and the field stronger, of the order 
of atomic fields. Then it will no longer make sense to talk of the motion of an atom in 
an external field. Instead we have an electron and a nucleus interacting and moving in 
a common external field. If, in addition, we fix the position of the nucleus we arrive 
at a situation similar to the single centre approximation, which is used in calculations 
of the P-odd energy shifts in chiral molecules [33|. There one considers single orbitals 
of electrons being exposed to the potential of one of the nuclei (single centre) and in 
addition to a chiral external electric field which summarizes the effects of the other 
nuclei. Again we see that chiral molecules and the situations we considered represent 
very different limiting cases in the above "Gedankenexperiment". 



We found that there are possibilities to enhance the P-odd effects in the energy levels 
dramatically, by many orders of magnitude. At a crossing of two levels of unstable 
states, i.e. complex energy levels in an electric field with nonzero handedness the P- 
violating energy shifts can become proportional to y/6i, i = 1,2. The comparison of 
these levels with the levels for the field configuration of opposite handedness shows 
a drastic difference: Under a reversal of the handedness the P-violating shift gets 
multiplied by a phase factor i. This means, that the shifts in the real and imaginary 
parts of the complex eigenenergies are exchanged. We emphasize that such energy 
shifts oc only appear for unstable states. 

Numerically we found the following results. For the case of a "hydrogen" atom 
with nuclear spin / = we obtained P-violating energy shifts proportional to \f8{ of 
AE/h ~ 10~ 2 Hz. For ^He" 1 " we found the same kind of shifts with AE/h ~ 1 Hz. 
Considering the real atom \H with nuclear spin / = 1/2, we found energy shifts 
proportional to \Jo~ 1,2 of the order of AE/h ~ 10 -5 Hz. These shifts are considerably 
smaller than the results of the calculation neglecting the nuclear spin of \B. because the 
narrow hfs-splitting restricted the reliability of our perturbative calculation to smaller 
electric field strengths. An improved calculation should overcome these restrictions 
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and could lead to much larger effects. In any case, our results show, that at the (would 
be) crossing points of the (P-conserving) complex energy levels, there remains a P- 
violating splitting and the change of this energy difference with the reversal of the 
field's handedness is enhanced greatly compared to the shifts away from the crossings. 

Finally we note that for the case of hydrogen discussed here all our energy shifts 
are much more sensitive to the parameter 82 than to 81 (cf. ( |1.10|) -( p..l2[ )). The reason 
for this is the same as given in Sects. 6.1 and 6.2 of Thus a measurement of our 
effects could lead to a determination of the nuclear spin-dependent weak charge Qff 
( p. . 1 1| ) which depends on the contribution of s-quarks to the proton spin. In this way 
parity violation in atoms could give information concerning the "spin crisis" of the 
nucleons. 
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Appendix 



A Representations of n=2 matrices and states 



In this appendix we discuss the non-hermit ian mass matrix A4 ( |1.39| ) and the dipole 



operator D as well as the matrix representations of the T and R transformations for 
atoms with nuclei of / = and 1/2 transforming like the a-particle and the nucleon, 
respectively. In addition we present the T- and R-transformation properties for the 
eigenvectors of the mass matrix. Here T is the time reversal operation and R denotes 
the reflection on the 1-3-plane. R can be considered as a product of a parity operation 
followed by a rotation by ir around the 2-axis: 







( Xi 




I:)- 








V ^3 


R = 







(A.l) 

All our conventions for the states etc. are as in fl5f . The matrices representing T and 



R are denoted by T and 71, respectively (cf. (3.33) and (3.58) of p"5fl .) 

The T- and R-transformation properties of Ai, of the quasiprojectors IP and the 
eigenvalues E (cf. Sect. |1.2| ; we drop the index INT for internal here) are given in Sects. 
3.3, 3.4, and Appendix B of flSfl . We have from T 

E(a,F 3 ,Se 3 ,8i,8 2 ) = E(a,-F 3 ,Ee 3 ,8i,8 2 ), (A.2) 

(l^W{a,F 3 ,Ee 3 ,5i,8 2 )T) T = E(a,-F 3 ,Ee 3 ,5i,8 2 ) (A.3) 

and from R: 

E(a,F 3 ,Ee 3 ,8i,8 2 ) = E(a, -F 3 , E e 3 , -8 X , -5 2 ), (A.4) 

E^E{a,F 3 ,Se 3 ,8 h 6 2 )n = E(a,-F 3 ,Ee 3 ,-8 h -8 2 ). (A.5) 



Since T = —1Z for / = and T = 1Z for / = 1/2, we get from (|A.2|) to (|A.5|) , combining 
the T- and R-transformations: 

E(a,F 3 ,Ee 3 ,8 x ,8 2 ) = E{a,F 3 ,Ee 3 ,-8 1 ,-8 2 ), (A.6) 
(P(a,F 3 ,^e 3 ,5 1 ,5 2 )) T = F(a,F3,5e3,-<Ji,-<y 2 ). (A.7) 

From the above relations for the quasiprojectors we turn now to the T- and R- 
transformation properties for the right and left eigenvectors of A4. Consider first the 
case of vanishing electric field, £ = 0. The corresponding eigenstates of JSdo are given 
up to order 8i, 2 in App. B, Tables B.IV and B.VIII of [|15| . For their T-transformation 
we have 

(P\2l,j,F,F 3 ,0e 3 ,8 1 ,8 2 )) T = (-l) F+F * +l (2i, j, F, -F 3 , 0e 3 , 5 1 , 8 2 \, 
(2l,j,F,F 3 ,0e 3 ,8 1 ,8 2 \ry = (—l) F+F3+l \2l,j, F, —F 3 , 0e 3 , Si, 8 2 ). (A.8) 
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Here I = 0,1,2,... or in spectroscopic notation S,P, ... denotes the orbital angular 
momentum quantum number and the hat indicates the perturbation by the P- violating 
Hamiltonian. Similarly we find for their R transformation 



]Z\2l,j,F z F 3 ,0e 3 ,5 1 ,5 2 ) 
(2l,j,F,F 3 ,0e 3 ,S 1 ,5 2 \n jf 



-1) F+F ^ l \2l,j, F, -F 3 ,0e 3 , -S 1} -8 2 ), 
-l) F+F *+ l (2l,j, F, -F 3 , 0e 3 , -5 lt -6 2 \ (A.9) 



By an argument similar to the one given below one can see that with a suitable nor- 
malization of the states eqs. ( |A.8| ) and ( |A.9| ) are valid not only up to order 5\ }2 but to 
all orders in <5i i2 . 

We will now define the eigenstates of Ai for £ ^ using the methods explained in 
App. C of 10. We set 

\a, F 3 , £e 3 , 8i, 8 2 ) 

= P(a, F 3 , £e 3 , 5 U 5 2 ) \a, F 3 , 0e 3 , 5 U 5 2 ) Ar 1/2 (a, F 3 , £e 3 , 5 1 , 5 2 ), 
(a,F 3 ,£e 3 ,8i,8 2 \ 

= (a,F 3 ,te 3 ,8 1 ,8 2 \W(a,F 3 ,£e 3 ,8 1 ,8 2 )Af- 1 / 2 (a,F 3 ,£e 3 ,8 1 ,8 2 ), (A.10) 



where 



F 3 , £e 3 , 8 h 8 2 ) = Tr [P(a, F 3 , £e 3} 8 1: 8 2 ) JP(a, F 3 , 0e 3 , 8 h 8 2 



(All) 



and a stands for (2l,j,F). For £ ^ these quantum numbers are not conserved and 
therefore indicate here the values for the corresponding state at £ — 0. From (|A.3 ) 
and (|A.8|) we now obtain immediately: 



Af(2l,j,F,F 3 ,£e 3 ,8 1 ,8 2 ) 
(P\2l,j,F,F 3 ,£e 3 ,8 l ,8 2 ) 



Af(2l,j,F,-F 3 ,£e 3 ,8 1 ,8 2 ), 
(-lf+^(2l,j, F, -F^£e 3 , 5 U 8 2 \, 

(-l) F+F3+l \2lJ, F, -F 3 , £e 3 , 5 U 8 2 ). (A.12) 



(2l,j,F,F 3 ,Se 3 ,5 1 ,5 2 \T 

Similarly we get from ( |A.5 ) and ( |A.9| ): 

N{2i,j,F,F 3 ,£e 3 ,8 u 8 2 ) = N{2l,j,F,-F 3 ,£e 3 ,-8 1 ,-8 2 ) 



£|2Z,j,i^F 3 ,£e3,<M 2 ) 
(2l,j,F,F 3 ,£e 3 ,5 1 ,S 2 \n^ 



-l) F+F * +l \2i,j, F, -F 3 ,£e 3 , -5,, -5 2 ), 
-l) F+F * +l (2l, j,F, -F 3 ,£e 3 , -5 U -6 2 \. (A.13) 



Up to now we have worked with the basis states (|1.41 ) which are simultaneous 
eigenstates of JSd (|1.39|) and F 3 . Alternatively, we may use the basis set of (left and 
right) eigenstates of Ai, which, for \F 3 \ ^ 0, are defined as follows: 

|a, \F 3 \,±,£e 3 ,S u S 2 ) 



i 

v/2 



oc . 



\F 3 \,£e 3 ,5 1 ,5 2 )±(-l) l ^- I n\a,\F 3 \,£e 3 ,-5 1 ,-5 2 ) 



(a, \F 3 \,±,£e 3 ,6i,5 2 \ 



i 

V2 



a. 



\F 3 \,£e 3 ,5 1 ,5 2 \±(-l) l ^- I (a,\F 3 \,£e 3 ,-5 1 ,-S 2 \]Z 



(A.14) 
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For F 3 = which occurs only for I — 1/2 we label the states as defined in ( |A.1U| ) with 
+ (— ) if F + l is even (odd). In ( |A.14j ) R denotes the reflection ( |A.1| ) on the 1-3-plane. 
For 5i = 5 2 = the states (|A.14|) are simultaneous eigenstates of A4 (|1.39|) and 71. 

As one can easily check, the definition of the relative phases between the two linearly 
combined _F 3 -eigenstates as well as the normalizations are chosen in such a way that 
the orthonormality relations 

(a, \F 3 \,r,£e 3 ,8i,8 2 \/3, \F 3 \, r' , £e 3 , 5 U 5 2 ) = 5 a/3 8\F 3 \\F!,\$rr> , (A.15) 

with r, r' e {+, — }, are fulfilled. The transformation under T, R, R°T of the states for 
the case / = is given by 



(l f \a, \F 3 \,±,Se 3 ,6 1 )Y = =r (a, \F 3 \, =f, £e 3 , S^, 



T 

r 3 |,±,£e 3 ,di|Tj T 

ti 

a,\F 3 \,±,£e 3 ,5 1 \U = T i (a, \F 3 \, ±,£e 3 , -5x1, (A.17) 



(a,|F 3 |,±,£e 3 AlZj = ±i|a,|F 3 |,T^e 3 ,(5i), (A.16) 
R: n^\a,\F 3 \,±,Se 3 ,6i) = ±i \a, \F 3 \, ±, £ e 3 , -tfA 



R°T: (|a,|F 3 |,±,£e 3 A)) T = (a, \F 3 \, T,£e 3 , 

((ajFslTi^es,^!) 7, = \a, \F 3 \, =f, £e 3 , —Si), (A.18) 

and for the case I = 1 /2 by 

T: (r^|a,|F3|,±,5e3,(yi,<J2)) T = ±{ot,\F 3 l±,£e 3 ,8 1 ,8 2 \, 

( K (a,\F 3 l±,£e 3 ,5 1 ,5 2 \T) T = ± \a, \F 3 \, ±, £e 3 , 8 U 5 2 ), (A.19) 

R: lV\a,\F 3 \,±,£e 3 ,5i,5 2 ) = ±\a,\F 3 \,±,£e 3 ,-8 1 ,-8 2 ), 

{a,\F 3 l±,£e 3 ,5 u 5 2 \ll = ± (a, |F 3 |7±, £e 3 , -S u -S 2 \, (A.20) 

R°T : (\a,\F 3 l±,£e 3 ,5 1 ,5 2 )) T = (a, |F 3 |7±, £e 3 , -8 h -5 2 \, 

((ajFglTi^es,^,^!) 7 = |a,|F 3 |,±,£e3,-<fi,-<y a ). (A.21) 

We note that for given a = (2/, j, F) and | ^3 1 7^ the corresponding + and — states 
are linearly independent and eigenstates of A4 to the same complex energy eigenvalue. 
For a = (2l,j,F) and F 3 = we have only the +(— ) state for I + F even (odd) as 
explained above, and there is no energy degeneracy. 



B The continuity matrix M c 

In this appendix we give details of the calculation which leads to the results of Section 
2.1 . We consider first, as in Sect. |2.1.1| , a neutral atom with an / = nucleus. 

Let us have a closer look at the eigenstates |(», /3) of the internal motion (cf. ( |2.13| )). 

In Sect. |1.2j we have recalled the properties of the right and left eigenvectors of the 



mass matrix (|1.39|) for an external electric field in 3-direction. Let us now consider 
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such electric field vectors of length £ (a) = (er = 1, 2), rotated into the 3-direction. 
As constructed in ( |1.41|) , we denote the right eigenvectors of J\A — D 3 £ (a) by 



|W, ft <$!,$,) := \a,F 3 ,£^e 3 ,8 u 8 2 ), (B.l) 

where (3 is a shorthand notation for (a, F 3 ), and we indicate explicitly the dependence 
on <5i 5 2- The corresponding left eigenvectors are 

By a proper rotation Ro we can bring the pair of vectors £ w and £ (2> into the 1-3-plane. 
By further suitable rotations around the 2-axis by angles 8 m or 6 (2> we can make either 
£ w or £ (2) to point along the positive 3-axis. With F,i (i = 1, 2, 3) the generators of 
rotations around the coordinate axes this is formally expressed as follows: 

exp [i6 m F 2 ] ■ Ro : £ m -> |£ (1) |e 3 , 

exp [i6 (2) F 2 ] ■ Ro : £ (2) -> |£ (2) |e 3 . (B.2) 



For the eigenvectors defined by ( [2.13 ) we can then set 

I W, (3, Sx, 8 2 ) = ©(Ro- 1 ) • exp \-i6^F 2 ] \{*)' , (3, 8 U 5 2 ), (a = 1, 2), (B.3) 

and correspondingly 

(H/3A,<* 2 | = ((*)', fc5 h 5 2 \ -exp^Fj -P(Ro), (a = 1,2). (B.4) 

Here 2? are the representation matrices of the rotation group in the space of the n = 2 
atomic states. From ( B.3|) , (p.4| ) we can write the matrices ( 2.15|) as: 



U^f(8 1 ,8 2 ) = (W, 7 , 6 U 5 2 \ exp [i - ^) F 2 ] |(r)', /3, 5 1; 5 2 ), 

(l<(7,r<2). (B.5) 



Thus we have expressed f/ (<TT) by the eigenstates (|B.1|) for electric fields pointing in 



3-direction and the relative angle 6 ,(CT) — 6 {r) between the fields £ {CT) and £ {T) . 

We will now show that ( p. 27] ) is equivalent to an equation having no linear terms in 
5\ Q2.28 ). In order to achieve this we choose the states defined in (|A.14[) of Appendix 



as a basis. In the notation \nl, F, |F 3 |, r, £e 3 , 8x) = \nlp, |F 3 |, r, £e 3 , 8x), (r = ±), 
for I = 0, j = F, as explained in (|A.10|) ff. of Appendix [A], we write the basis of right 
eigenvectors as: 



251/2, 


1 

2 


r, 


\£ {a) 


|e 3 ,5i 


I2A/2, 


1 

2 


r, 




|e 3 , 81 


I2P3/2, 


1 

2 


T 


\£ {a) 


e 3 , <$i 


I2P3/2, 


3 
2 


r 
1 1 


\£ {a) 


|e 3 ,5i 



±,(7=1,2) (B.6) 
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and the analogous set for the left eigenvectors. We now write a, r, Si), a = 1, 2, 3, 4, 
as shorthand for the states ( |B.6|) and define the following 4x4 matrices (cf. ( |B.5|) ): 



(IZ (CT ' T)±± (5i)) Q/3 := ((«r)' l a,± l tf 1 |e i ^^|(r/ )i 9 ) ±, < 5 l ) ) 
(U^(5i)) a p := ((<^=Mi|e i(CT - T) ^|(</3,TA), 

where a, r G {1, 2} and 

e = # (2) - e (1) 



(B.7) 
(B.8) 

(B.9) 



is the relative angle between the fields £ (2) and £ m . From the transformation properties 
of F 2 and of the states (|B.6|) under R (cf. ( |A.18|) ) we find 



(B.10) 



(B.ll) 



W"' T)rs (Si) =r ■ s ■U^ T)rs {-5 1 ) 

and from the R°T transformation (|A.18 ) we get 

[U (a - T)rs (5i)] T = f/ (T " T)( - 3)( - r) (-5 1 ). 

Therefore, as a consequence of the T-transformation, the [/-matrices obey 

\U<" r *'(6i)] T = r-s- £/ (T ' CT)( - 3)( - r) (5i). (B.12) 

From ( |B~T0D we see that U (a ' T)±± in (|R7]) is even, U (a < T)±T in ( gj) is odd under ^ -> 

Let us furthermore define 4x4 matrices S {a) , C^' as in (|2~23|) . Note that the values 
Qa defined in ( P-17| ) are the same for the + and — states. Putting everything together 
we find for M C (E, Si), in the basis of the ± states and with a suitable arrangement of 
rows and columns: 

M e (E,6 1 )=M' e (E,6 1 ).MZ(E,6 1 ), 



where 



ML 














[|(1,2)++ 



















£/(l,2)- + 







(B.13) 



(B.14) 



M':(E,Si) 



( s w 











\ 



















s w 







V o 








_ C (2) 


) 



From (|B.13[) we get 



(B.15) 



(B.16) 



det M C (E, Si) = det M' C {E, 8 X ) ■ det M%(E, Si). 



(B.17) 
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For general fields £ (1> , £ (2> we have det M" ^ and therefore det M C (E, 5\) = implies 
det M' C (E, 5i) = and vice versa. Now we know that in M' c (P-14Q the submatrices 
Lf (<T,T)±± are independent of, and the U l ' 7 ' T)± " F are linear in Si, neglecting terms of order 
Sf. Thus we have 

det M' C (E, 5 t ) = det ^ ^ c &1 £ J = det{CAC- x B - 5fCD), (B.18) 

where A, ...,D are square matrices depending on E but not on 5\. Therefore we can 
write the equation determining E as in ( ]2.28|) . 



Now we will discuss the properties of the matrix M' c ( |B.14| ) and give the proofs of 
(pT29|) (1273T1) . 



After rearranging rows and columns of M' c and multiplying some columns with 
'— 1), we find 

det M' C (E, 8x) = det M? sym {E, 6J, (B.19) 

r(2,i)++ t (2) 



( 



\ 





juy-A- 
-Tg*> 





-a (l) 



U {2 ' 1)++ Tg m \ 

Q<X> f/(l,2)++ 

u {1 ' 2) ~ + 

I7(2,i>-+ o / 



(B.20) 



Using now the T relations (|B.11| ) and the definitions ([2.231) , ( |B.15|) of the matrices 
Tg {,7) , Ct (rT \ it follows immediately that the above matrix is antisymmetric. Note that 
Jfasym ig a 16 x 16-matrix. 



A basic theorem of linear algebra — that the determinant of an antisymmetric 
matrix M with an even number of rows is a square — may be applied now. In order to 
prove this theorem, one writes the matrix M = —M T , assuming M ^ (after possibly 
necessary rearrangements of rows and columns) as 



M 



A B 
-B T C 



with A :- 



a 
-a 



a ^ 0, 



(B.21) 



being a 2x2 matrix and C = —C T antisymmetric. Then, with 



detM = det{( B 4, Ma/ (J ~ A ;' B 



det 



A 

C + B T A- 1 B 



(B.22) 

it follows that det M = a 2 det M', with M' = C + B T A~ X B being antisymmetric. The 
theorem now follows by mathematical induction over the number of rows of M. 



Applying this theorem to M^ sym ( |B.20|) and using ( [2.28] ) we find: 
det M c {E,8 l ) = det Mr m (E,5 1 ) = [g(E,5 1 )] 2 = f {E) + 6 2 J 2 {E) + 0{5f), (B.23) 

g(E, 5 t ) = g (E) + 8 2 g 2 {E) + 0(«Jf). (B.24) 
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From ( |B.23p we see that the energy eigenvalues E are at least twofold degenerate for 
all fields £ (1) ,£ (2) up to order 6\ (Terms of higher order in 6\ have been neglected in 
previous stages of the calculations). We know that for = £ (2) there is at most a 
twofold degeneracy, thus, apart from possible special values of £ (1) , £ (2) the degeneracy 
will not be bigger than twofold. It follows that for generic values of £ {2) go(E) will 
have only non-degenerate zeroes. The energy eigenvalues are now obtained from 



g(E,5 1 







(B.25) 



and using the implicit function theorem we find from (|B.24|) that the solutions E 
E{5%) of ( |B.25| ) have no terms linear in 8\\ 



dE(5 x 



0. 



(B.26) 



5i=0 



This completes the proof of the assertions made in Sect. |2.1.1| for the case of nuclear 
spin 1 = 0. 



Finally we have to deal with the case I = 1/2 and to prove that ( [2.33P is equivalent 
to ( ggg) (cf. Sect. |Xg ). 

Using the ± eigenstates as defined in ( |A.14| ) and the subsequent text of Appendix 
0we introduce matrices U {a ' T)ra (6 1 , 6 2 ) (r, s G {+, -}), as in (|B77D , (|B~8|) , for which the 
transformation properties under R, R°T and T analogous to ( |B.1(J| )-( |B.12| ) now read: 



R 

R°T 
T 



U 



(cr,r)r s 



(61,62 



r ■ s ■ U 



(cr,r)r s 1 



-61, -6 2 



[U^ )rs (6 1 ,5 2 )] T = U^i-St, -8 2 ), 



[U^ a (8 1 , 6 2 )f = r-s- U^ )sr (6 u 6 2 ). 



(B.27) 
(B.28) 
(B.29) 



From ( |B.27| ) we see again that U {tT ' T)rs is symmetric (antisymmetric) under (61,62) 
(—81, —62) for r = s (r s). The matrix M c factorizes as follows: 



with 



M C (E, 6x, 62) = M' C (E, 6 1 , 62) ■ M^(E, S u 6 2 ), 



r(2,i)H 



M' 



a (1)+ 





Ml 



( S w 



V 












£/(2,lH 



f/(2,D- 

a (1) - 




s w - 






f/(1.2)4 

~Tf y 





_ Ci 2 h 



I 



\ 



(B.30) 



(B.31) 



(B.32) 



J^(2)r _ <jf(2)r/£f(2)r\-l 

a (1)r = c mr (s mr )-\ 

r e {+,-}. 



(B.33) 
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Here S (2)r , C (1)r are the diagonal matrices as defined in ( |2.23| ) for the r = +, — states. 
Now the important difference to the 1 = case is that for I = 1/2 the r = + and 
r = — sectors are not simply related since the F3 = states with r = ± are not energy 
degenerate (cf. Appendix A). Thus the corresponding values Qp in ( f2.23| ) are different 
and also S {2)+ ^ S {2) ~,C W+ ^ C w ~. 

We can draw the following conclusions for the eigenenergies E obtained from 



det M' C (E,5 1 ,8 2 ) = 0. 
For 5\ = 62 = the matrix M' c is block-diagonal since 

C/ (2 ' 1)±T (0,0) = U {ia)±T (0,0) =0. 



(B.34) 



(B.35) 



But for general fields £ (1) ,£ (2) this does not lead to an energy degeneracy, since the 
++ and — blocks of M' c have different determinants due to the non-degeneracy of the 
F 3 = plus and minus states. For 5±, 5 2 unequal zero, we see from ( P-31| ) that we can 
apply the same arguments as in (|B.1£|) and 



in order to obtain (|2.34). 



C A modified version of Bloch's method 
for degenerate perturbation theory 

In order to calculate eigenvalues and eigenvectors of a time independent operator O in 
a perturbative approach, where O is a sum 

d = d + xd 1 (c.i) 

of an operator Oo, whose eigenvalues and (orthonormalized, complete set of) eigen- 
vectors are known, and a perturbation XOi with a small parameter A, one uses time 
independent perturbation theory to calculate perturbative corrections to these eigen- 
values and eigenvectors. If (some of) the unperturbed eigenvectors of O have the same 
eigenvalues, which e.g. is the case in most of our problems, one has to use degener- 
ate perturbation theory (cf. e.g. For this case of degeneracies, several methods 
have been given to write down the explicit perturbative expansion of the eigenvalues 
of 6. In Chap. § we have used both Kato's (cf. (13371) -(ft^IQ) and Bloch's (cf. (P^)- 
( p.31| )) method. The advantage of Bloch's method is that by modifying H (cf. Q3.27] )) 
one obtains an eigenvalue problem ( |3.31| ), where the sought eigenvalue E is multiplied 
only by a unit operator .0 This eventually yields simpler expressions (compare ( |3.28|) ff. 



with (|3.17|) ff.). However, if H and K are hermitian H in general is not because of the 



asymmetric multiplication with K in ( p. 27 



HK = H. (C.2) 



3 Since the discussion here is general we drop the indeces a or (a/3) on the operators H, K, as 
well as on the eigenvalues E, which in the main text are used to indicate the subspace of unperturbed 
states in which the diagonalization takes place. 
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In order to obtain both, an eigenvalue problem of the form ( |3.31| ) and hermiticity, we 
introduce a modified version of Bloch's method (cf. (|3.78|) f.), which we describe in 
more detail in this appendix. 

The method is similar to the well-known treatment of small oscillations in classical 
mechanics, where the equations of motion, after a Fourier transformation, lead to the 
eigenvalue problem 

(if-wjT)-a* = 0. (C.3) 

Here V = (Vy) and T = (T^) are the second order coefficients of the expansion of 
the potential and the kinetic energy about the minimum of the potential w.r.t. small 
displacements of the generalized coordinates. The elements of the eigenvectors give 
the Fourier coefficients of the tranformed solutions of the equations of motion with 
eigenfrequencies Uj. (cf. e.g. |3"1"H). 

If T is hermitian and positive, it possesses a hermitian and positive inverse T~ l and 
square root T 1 / 2 with (T 1//2 ) 2 = T. Then (|C.3|) is equivalent to 

(T-VVr-i/2 - 4) ■ a' fe = 0, (C.4) 

with a' k = T 1//2 afc. Therefore it suffices to find a set of eigenvectors a' k of U := 

We extend this method to our case of the (in general) non-hermitian matrices H 
and K (the analogues of V and T respectively) and define a matrix H by 

K^HK 1 ' 2 = H, 
ffP (0) = H, 

P m H = H. (C.5) 
Our eigenvalue problem ( |3.21| ) which is similar to ( |C.3| ) is then equivalent to ( |3.79| ): 

(H-E)\ ) = 0, (C.6) 

where | ) G TZ (0) . 

From the expansions of K and H, Q3.17 ) ff., we find the expansion of i^~ 1//2 , the 
inverse of K 1 ' 2 in the space 1Z {0) . This inverse exists for sufficiently small £', since K 
is then arbitrarily close to P (0) which acts as the unit matrix over TZ m (cf. ( |3.17 )). Up 
to fourth order in £' we get: 

£-1/2 = pw + I(o,2,0) 

+ \ [(0, 1, 2, 0) + (0, 2, 1, 0) + (0, 3, 0, 0) + (0, 0, 3, 0)] 

+ §E (4) (Mi, ^3,0) -§(0,2,0,2,0) 

+ 0{£' 5 ). (C.7) 
From (|3.17| ) f. and ( |C7| ) we calculate the expansion of H as 

H = k~ 1/2 HK- l/2 



£(°)pt°) + £-1/2 ^B ,(n) j k- 1 ' 2 

oo 

£(o)p(o) _|_ ^ C {n \ (CA 



n=l 
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with the terms of order n 

C (3) : 
Q(4) : 



1,2,3,4: 



(0,0), 
(0,1,0), 

(0,1,1,0) + | [(0,2, 0,0) + (0,0, 2,0)], 
(0,1,1,1,0) 
+ \ [(0, 0, 1, 2, 0) + (0, 1, 2, 0, 0) + (0, 2, 0, 1, 0) 
+ (0,0, 2, 1,0) + (0,2, 1,0,0) + (0,1, 0,2,0) 
+ (0,0, 0,3,0) + (0,3, 0,0,0)". 



(C.9) 



Finally we will discuss the T-invariance properties of Kato's eigenvalue problem 
and show that the matrices H aSe , Hss and H & sHg & obtained by our modified version 
of Bloch's method (cf. ( |3.80| ), ( |3.83| ) and Section |3.1| ) are all multiples of the 2x2 unit 
matrix. 

read: 



Let us consider first Kato's operator K. The matrix elements of Kt & p\ 

(^(^))(7',m'),(7,m) = if, m^K^g^f, 



(CIO) 



where 7', 7 G {a, /3};m' ,m G {+1, —1}- For the diagonal block matrices with 7' = 7, 
the same T-invariance arguments as in ( |3.23| )-( |3.25| ) apply, which yield that they are 
multiples of the unit matrices in IZ^.IZ^] cf. ( |3.26|) . As far as the off-diagonal matrices 
7' 7^ 7 are concerned, it follows from ( |A.12j ) that 

(f^n'\K m \^m) = (-l)^'^+^+™'+HI^I ( 7~ m |^ ((5 - ) | f) _ m0 . (C .n) 

Analogous relations apply for the matrix elements of Hr a j3)' 

We can now explicitly check that e.g. H a ^Hp a is a multiple of the 2x2 unit matrix. 
From (|C.11|) and the analogous relation for Hi^p)-, together with ( |3.78|) , we get: 



(a,m\H a pHp a 



a, m 



(a, m\ H {ccp) I A rn')(P, m'\H^\a, m") 

m'=±l 

V" (_iy{F+F'+m'\F 3 \)+(m+m")\F 3 \ 
m'=±l 



(C.12) 



x a, 



-m 



»\H m \P,-m')(fi, -m'\H { 



a/3) I 



Since 2(F' + m'\F 3 \) is even, it follows that IF + (m + m")]^ is even/odd for m" = 
m / m" = —m, respectively. Because we sum over m' in ( |C.12| ) this completes the 
proof. To show that H a ^(H^ — E)H^ a is also a multiple of the unit matrix, we still 
have to use the analogue of ( |3.25| ) for Hm- 



D General proof of the existence of y^i-resonances 

In this appendix we give a proof that for sufficiently small £' there exist "resonance" 
values for the vectors r7 (,T) where energy shifts proportional to v^i occur (cf. (|3.95|) - 
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Let us again first consider the case of a "hydrogen" atom with nuclear spin 1 = 0. 
We expand a and b ( |3.88|) in powers of £' up to second order and c, which has the 
dimension of a 2 ,b 2 , up to third order, where <5i-linear contributions of lowest order in 
£' show up. This gives 



a 
b 

c 



a^ + a^£' + a^£' 2 + 0(£' 3 ), 
b^ + b^£' + b^£' 2 + 0(6' 3 ), 



(D.l) 



All coefficients al n) , 6- n) , c-™ ) , i = 0,1; n = 0,1,2,..., depend still on the t7 (ct) and of 



course on the unperturbed configuration parameters (|3.61|) , (|3.68|) , ( |3.69| ), but not on 
£'. From the definition ( [3. 78j ) of Hr & p\ together with (|3.76|) and (|3.77|) it follows, that 
(cf. (FB5D ) 



°0 



77.(0) 



0. 



The series for c in (|DlD impl les 



-(3) 



+ ^ci- 



cs) 



£' 2 + 0{£ 



I3\ 



(D.2) 



(D.3) 



As one can see, the third order term of c is needed in order to obtain a consistent 
expansion of ±i(a — b) + l^fc up to second order in £'. Since we = ( (ft.66|) , 
( p.67| )), the lowest order contribution to the radicand in ( ft.87Q will be oc £' 2 . Therefore 
we can write the l.h.s. of ( p.91|) , ( |3.92| ) as follows: 

±i(a -b ) + 2^ =: £' ■ F±({r,^}, £') 



±i( a w _ #)) + 2 Jc^ + 0(6') 



£'. 



(D.4) 



Suppose now a set of perturbation fields ( |3.97| ) can be found satisfying Q3.9H ) or ( |3.92| 



at order £' i.e. 



F + ({77< CT >},£' = 0) = for 17 W = f&U (D.5) 
F„({ V ^},£' = 0) = for rj^ = fj^ s> _, (D.6) 

where also ( p.93|) holds. Considering e.g. ( p.5|) we can then apply the implicit function 
theorem in order to show that a solution to ( |3.91| ), i.e. 

F + ({^},£') = (D.7) 

for some values T7 (CT) = »7fe S -f exists to all orders in £' for sufficiently small £'. 

Now the condition that the implicit function theorem may be used is that (a) F + 
is a continous function and that (b) the functional determinant 







drii 













7^0 



(D.8) 



,res ,+ f^ij 
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with respect to some pair of coefficients 77^ ^ r/j T> does not vanish at the zero of F + (£' = 
0) which is determined by (p.5|) . Similar statements hold for F_. We emphasize that we 
can in this way give a general existence proof of P- violating energy shifts proportional 
to y/S\ using only the lowest relevant orders of perturbation theory. 

For the case of our "hydrogen" atom in an unperturbed configuration described 
by ( |3.61| ), ( |3.68| ), ( pj.69| ), the following perturbation field configurations satisfy ( p.5[ ), 

& 



, (D.9) 





( ^ \ 








f ^ \ 


'/res,± 






I " ) 


j '/res,± 






V ±77.08 fj J 




V ^480.58r) / 




\ ±77.30 fj J 



where 



^=V(C-C) 2 + ^ 1) - J |f) 2 (D.10) 

In general, there will be a continuous set of solutions of (|D.5| ), (|D.6|) . In order to 
simplify the calculation, we have chosen our solutions such, that they also satisfy 

&(4 1} - #>) = 0. (D.ll) 

The functional determinant (p.8| ) with respect to 773 — rj^ and 773 + r/ 3 3) evaluates to 



(3) (1) 1 (3)\ 



T 2.63 • 10 _5 (MHz h V _1 cm) 



(D.12) 



This completes the proof, that to any order in £', there exist non-trivial solutions of 
( p.91| ) and ( j3.92j) . In general one expects, that for such a solution also ( |3.93j ) holds. 
We will now give explicit solutions of (|3.91|) , (|3.92|) to second order in £', which yield 



a P- violating energy shift proportional to V"i (|3.94|) . 

We expand the two P-even factors of the radicand in ( |3.87|) , i.e. ±z(ao — 60) + 2y/c~o 
up to second order in £' (cf. ( |U) . 1| ) , ( p.3| )): 

± t(a - b ) + 2^F =: £' • ({17^}) + F^ii^})? + 0(£' 2 ) 



±t{a ( 1) 



+ 2\M 



£' 



+ 



±t(a ( 2) 



£ 



l2 



+ 0{£' 3 ). (D.13) 

Here F^\{r} (a) }) = F±({r) (a) } , £' = 0) and Fj. l) ({i7 w }) are the zero and first order 
coefficients of the expansion of F± ( P-4j ) in terms of £'. Now, for £' 7^ and to second 
order in C |379T|) and ( ggg ) imply 



F { ±\W a) }) + F£\{q w }) ■ £' = 0. (D.14) 

Again, for our "hydrogen" atom in an unperturbed configuration described by ( |3.61|) , 
( p.68| ), ( |3.69D , the perturbation field configuration Q3.62| ) satisfies ( p.!4[ ) with either: 

/ 


V -702.30 



'Ires, 



1 



117.00 



'Ires, 



„(3) 
'Ires, 




;d.is) 
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for 



or 



„Ci) 

'Ires,— 





„(3) 
'Ires,— 




(D.16) 



(D.17) 



for 

£' = 3V/cm. (D.18) 

Note that the perturbation fields ( p.!6| ), ( |D.18| ) have been chosen in a way that they 
fulfill the condition \£'\ < £' max , where £' max is obtained by an argument as in ([3.64]) 
adapted for our case here. This gives an estimate of the maximum for which the 
perturbation series is reliable. £' max is found to be 



£' 



9V/cm, 



£' 



12V/cm 



(D.19) 



for the "+"- and "—"-solutions ( p,15[ ) and ( |L).17| ) respectively. Within these limits, 
we looked for solutions which for simplicity we required to satisfy: 



ft (a« - 6" + 5'[4 2) - &8 ]) = 0. 



(D.20) 



This led to further restrictions on £', and the values given in ( |D.16| ), ( |D.18j ) represent 
approximately the maximal ones compatible with (|D.20| ). Then for the "+" -solution 
( p,15| ) £' turned out to be an order of magnitude smaller than for the "—"-solution 

(ECUD- 



To finish we turn to the case of an ion with 1 = 0. For the case of a ^He 4 " ion in a 
box described by ( J3.61[ ) , ( |3.109| ), ( |3.11U| ), the following perturbation field configurations 
satisfy O, (TO: 



T7 (1) + 
'Ires,± 



- (2) 
"Hres,± 



r? (3) + 

'Ires,± 



( 



\ 



\ 0.02*71 =f (62.02^0.11; 






V1 + V2 J 



y -0.06 fjx ±445.77^ + 7^ j 



\ 0.0277! =f (62.02 ± 0.11) yjfft +f% J 



(D.21) 



where 



_ „(3) 

71 <h > 



(3) 



(D.22) 
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In general, as for q = 0, there will be continuous sets of solutions of ( |D.5[ ), ( P-6j ). In 
order to simplify the calculation, we have again chosen our solutions such, that they 
also satisfy (p.ll[) . 

The functional determinants ( p.8| ) with respect to r]^ — rjf- 1 and rj^ + ?73 3> evaluate to 



^4.70 • 10 -8 MHz h V _1 cm 



(D.23) 



This completes the proof for ^He , that to any order in £', there exist non-trivial 
solutions of ( |3.91| ), ( |3.92| ) (cf. the discussion concerning the implicit function theorem 
in |3J]). For such a solution one in general expects, that (|3.93|) holds as well. As in the 
previous section, we will now give explicit solutions of ( |3.91| ), ( |3.92| ) to second order in 
£', which yield P-violating energy shifts proportional to \f5[ ( |3.94 ). We proceed as in 



( p,13| ) f. The following perturbation field configurations satisfy ( p.l4| ): 









( 








„w _ 

'lres,+ 


I 1 ) 


j 'lres,+ 




1 ) 


„(3) _ 

) 'lres,+ 


[ ' ) 




^ 94.45 / 




\ 


-653.30 / 




V 94.62 / 



(D.24) 



for 



and 



£' = 9.6V/cm 



;D.25) 



'Ires, 




'Ires, 






V 356.97 



„(3) 
'Ires, 




(D.26) 



for 



£' = 96 V/cm. 



(D.27) 



As in the q = case the perturbation fields ( p.25|) , ( p.27|) satisfy the condition \£'\ < 
£' max , where £' max is given by ( |3.64| ) and is an estimate of the maximum for which the 
perturbation series can be used. £' max is found to be 



£L 



192 V/cm, 



£L 



320 V/cm 



'D.281 



for the "+"- and "—"-solutions ( P-24|) and (p.26|) respectively. Within these limits, 
( p.25|) and (p.27|) are chosen such that solutions exist, which for simplicity we again 
chose to be restricted to those which satisfy (p.20|). 
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Table Captions 



Table I. The P-conserving as well as the P-violating perturbative contributions to the 
eigenenergies Ea,± (|3.128| ) ff. up to 2nd order perturbation theory for different lev- 
els a = (a, I-F3 1 ; n = (1, 1, 1)) of }H at different electric fields £ for the corresponding 
maximally allowed perturbation field strengths S' max . 

Figure Captions 

Fig. 1. Sketch of the electric field configuration in a box divided into three segments (l)-(3). 
The electric field is homogeneous within every segment and varies suddenly along the 
1-direction. 

Fig. 2 a,b. The energy levels = E m {a, |F 3 |,£e 3 ;n) (|3.10|) for a "hydrogen" atom with 
/ = in an external electric field <?e 3 and n = (m, n 3 ). The real parts for a = 2S\/2 
and 1P\j2 and n 2 = n 3 = 1 are shown in a, the decay widths (twice the negative of the 
imaginary parts in units of s _1 ) for a = 2S , 1 / 2 , 2P1/2 an d 2P 3 / 2 in b. In our calculation 
the imaginary parts are independent of n. 

Fig. 3. The real and imaginary parts of the P-violating energy- difference AE^ -£ 3 /(r^' 3 ^) 
(cf. ( |3.63D ) as a function of £ for the state a = (a, |F 3 |; n) = (2S' 1/ / 2 , 1/2; (1, 1, 1)) within 
a range of £, where extrema due to a level crossing occur (cf. Fig. 2a). 

Fig. 4. The field strength value £' max characterizing the limit of applicability of our pertur- 
bative treatment (cf. fl3.64|) ) as a function of the field strength £ for the same level a 
as in Fig. 3. 

Fig. 5 a,b. The P-violating energy shift AE^ ( |3.63j ) for £' = £' max - The real part is plotted 
in a, twice the imaginary part in b. 

Fig. 6 a,b. The decadic logarithms of the moduli of (a) the real and (b) twice the imaginary 
parts of the P-violating energy-difference E^s + ({£' (ct) }) — E & p + ({£ H (<T) }) ( |3.99| ) for our 
"hydrogen" vs. the deviation 773 — r]^l es _ of ry 3 2) from the resonance value r]^l es _ = 
574.9 for which the P-even splitting of levels is removed, revealing a P-violating splitting 
of the order of ■\f5\- 

Fig. 7. The lowest energy levels of the cm. motion of a 2 He + ion moving in one dimension 
in a linearily rising electric field potential of a field strength £ within an interval of 
length A = 1 nm. These are the eigenenergies of the Airy states (cf . ( [2.43D ) . For £ = 
they show the quadratic behaviour of the eigenenergies of sinus states (cf. ( |2.11D ), 

Fig. 8. The field strength value £' max characterizing the limit of applicability of our per- 
turbative treatment of 2 He + as a function of the field strength £ for the level 
a ={2S 1/2 , 1/2; (1,1,1)). 

Fig. 9 a,b. The P-violating energy shift AE^ in 2^- e+ f° r £' = £'max f° r the same level a 
as in Fig. 8. The real part is plotted in a., the imaginary part in b. 

Fig. 10 a,b. The decadic logarithms of the moduli of (a) the real and (b) twice the imaginary 
parts of the P-violating energy-difference E 5l ^ + ({£' (a ''}) — E^ +^E^}) ( |3. 113| ) for 
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2 He + vs. the deviation r/ 3 3) — rjs jres - of 7? 3 3) from the resonance value rj^ res _ = 525.4 
for which the P-even splitting of levels is removed in 2nd order perturbation theory, 
revealing a P-violating splitting of the order of \f5\. 

Fig. 11 a,b. The decadic logarithm of the moduli squared of (a) the real and (b) twice 
the imaginary parts of the P-violating energy-difference Ea,±({£ /( - a) }) — Ea,±({£^}) 
( 3.146| ) for }H vs. the deviation <? 2 (3) — £ 2 (3> e <, of £ 2 (3) from the resonance value £ 2 ( rL = 
—1.64 V/cm for which the P-even splitting of levels is removed in 2nd order perturbation 
theory, revealing a P-violating splitting of the order of y/fa. 
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